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SECTION  1 
Introduction 


In  recent  years,  there  has  been  an  increased  interest  in  the  electromagnetic  scattering 
from  forests.  Current  measurements  of  the  scattering  from  trees  cover  almost  the 
entire  electromagnetic  spectrum  from  low  frequency  radio  waves  to  light  [1,  2].  One 
frequency  band  of  interest  is  in  the  low  microwave  range  near  L-Band  spediically 
between  0.3-1.3GHc.  This  frequency  range  is  of  interest,  because  the  electromagnetic 
wave  is  of  long  enough  wavelength  to  penetrate  through  the  outer  leaf  canopy  of  the 
forest  yet  is  small  enough  to  provide  information  about  what  is  below  the  leaves. 
These  characteristics  make  these  wavelengths  useful  in  remote  sensing  applications. 
Previously,  most  of  the  measurements  and  modeling  of  the  scattering  from  forests 
at  L-Band  has  concentrated  on  narrow  band  radar  systems.  These  systems  have  low 
resolution  and  are  unable  to  isolate  the  scattering  from  individual  trees.  Therefore 
when  modeling  is  needed,  the  scattering  behaviour  from  an  individual  tree  is  not 
imp<ntant. 

For  nltrawide  bandwidth  radar  systems,  the  scattering  from  a  singfle  tree  is  im¬ 
portant  since  the  individual  trees  are  identifiable.  Therefore,  the  electromagnetic 
scattering  from  an  individual  tree  needs  to  be  examined.  To  study  the  scattering 
from  a  tree,  a  modified  form  of  the  Geometrical  Theory  of  Difiiraction  (GTD)  for 
dielectric  materials  is  used  for  the  analysis.  For  background,  an  introduction  to 
GTD  is  presented  in  Chapter  2  to  provide  the  basic  information  used  in  analyzing 
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the  scattering  from  a  tree.  Chapter  3  provides  background  into  current  methods  and 
parametos  used  in  modeling  trees,  wood,  and  soil. 

Before  the  geometry  of  a  tree  is  analyzed,  a  simpler  geometry  is  examined.  To 
do  this,  measurements  of  the  radar  cross  section  (RCS)  of  tree  logs  are  presented  in 
Chapter  4.  The  measurements  of  the  logs  help  to  conJirm  the  methods  and  param¬ 
eters  used  in  modeling  wood  that  are  presented  in  Chapter  3  since  the  cylindrical 
geometry  of  the  logs  is  simpler  to  analyze  than  the  complete  geometry  of  a  tree. 

After  the  wood  modeling  is  confirmed  through  log  measurements,  the  scattering 
from  various  parts  of  a  tree  are  examined.  Chapter  5  examines  the  scattered  fields 
in  the  vicinity  of  a  tree  trunk.  For  this  analysis,  the  trunk  is  modeled  based  on  the 
results  of  Chapters  3  and  4.  In  Chapter  6  a  model  for  the  scattering  from  a  tree  that 
includes  interactions  between  the  ground  and  the  trunk  are  examined.  In  Chapter  7 
the  scattering  from  tree  branches  is  examined.  By  separating  the  scattering  from 
the  trunk  and  ground  from  the  branch  scattering,  the  contributions  and  importance 
of  each  scattering  mechanism  is  determined. 
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SECTION  2 


Geometrical  Theory  of  Diffraction 

1  Introduction 

This  chapter  provides  a  brief  description  of  the  modified  Geometrical  Theory  of 
Difiaction  (GTD).  The  £ar-sone  diffracted  fidd  and  diffraction  coeffidents  for  per¬ 
fectly  conducting  and  didectric  wedges  are  presented.  Also,  corrections  to  GTD  for 
curved  disks  are  shown.  The  term  GTD  is  used  for  convenience  but  is  meant  to 
incorporate  the  Uniform  Geometrical  Theory  of  Diffraction  (UTD)  and  the  modified 
form  of  GTD  for  didectric  media  throughout.  Details  on  GTD  can  be  found  in 
Ite£etences[3, 4,  5].  The  time  dependence  e^  is  used  and  is  suppressed  throughout. 

2  Theory 

The  Geometrical  Theory  of  Diffraction  (GTD)  was  devdoped  in  great  detail  by 
KeDer  and  can  be  found  in  [6].  The  general  solution  to  a  problem  has  the  form: 

E‘(P)  =  M\P)  +  i'iP)  +  E^{P)  (1) 

where  ^(P)  is  the  inddent  fidd  at  the  observation  pdnt  P  and  M'^P)  is  the  reflected 
fidd  at  P.  The  far  none  diffracted  fidd  .^(P)  at  point  P  is  of  the  form: 

=  (2) 

O 
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where  7  ii  a  dyadic  dif&action  coefficient  pven  by: 


(3) 


and  .£*((?«)  is  the  incident  field  at  the  point  of  diffiaction  Qe.  The  amplitude  spread¬ 
ing  &ctor  is  defined  as  the  square  root  term  in  (2).  For  the  perfectly  conducting 
case,  the  diffraction  coefficients  D«  and  Da  corresponding  to  the  soft  and  hard  po¬ 
larisations;  respectively,  are: 


The  caustic  of  the  diffracted  ray  is  and  is  given  by: 

1  _  1  I  n,»(i'  -fs) 

ff^  p,  pgiaxi^fio 


(5) 


where  pe  is  the  radius  of  curvature  of  the  inddent  wave  in  the  edge  fixed  plane  and 
Pg  is  the  radius  of  curvature  of  the  edge  at  the  point  of  difhaction.  The  geometry 
for  these  diffracted  rays  is  shown  in  Figure  1. 

The  above  di&action  coeffident  in  Equation  (4)  must  be  modified  for  the  didec- 
tiic  cases  that  are  encountered  in  analyzing  trees.  It  is  shown  in  References  [7,  8,  9] 
that  the  modified  di&action  coeffident  for  an  dectrically  thin  didectric  is: 


2nV^irfcsin)9o 

{ |(i  -  “*  (--I--)  +  (1  -  + 


where  and  are  the  reflection  coeffidents  for  the  soft  and  hard  polarizations 
for  the  V  and  the  ‘n’  faces,  respectivdy,  and  T***  and  T'**  are  the  transmission 
coeffidents  for  the  soft  and  hard  polarizations  for  the  *o’  and  ‘n*  faces,  respectivdy. 
Calculation  of  the  reflection  and  transmisrion  coeffidents  are  pven  ir  \ppendix  A  or 
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Figure  1:  Canooical  geometry  for  the  di&action  by  a  wedge. 


in  Reference  [10].  The  above  equation  can  be  reduced  further  if  the  assumption  that 
the  dielectiic  material  is  lossy  so  that  the  transmission  coefficient  becomes  small. 
>^th  this  assumption,  the  diffraction  coeffident  reduces  to: 

[JJ.*  cot  +!!;•*  cot  (7) 

This  equation  can  then  be  used  to  calculate  the  diffraction  from  lossy  dielectric 
media. 

3  Chiang  and  Marhefka  Rim  Caustic  Correction 

The  above  formulas  for  the  diffracted  field  don’t  provide  the  correct  diffiacted  field 
for  the  drculai  dish  on  the  ends  of  a  cylinder.  To  produce  the  proper  diffracted 
fields  from  the  end  face  of  the  cylinder,  a  rim  caustic  correction  factor  needs  to  be 
induded  m  the  diffraction  coeffident  to  produce  the  proper  fidds  from  the  cylinder 
ends.  An  appropriate  caustic  correction  for  the  cylinders  encountered  in  anal]rsing 
trees  is  |pven  by  Chiang  and  MarLefka  in  [11].  The  new  diffraction  coeffidents  for 
the  perfect  electrical  conducting  (PEC)  case  are  pven  by: 


where:  _ 

r.(t7)  =  (») 

and: 

A(P)  =  x[^-fl,(l7)  +  M(I7)].  (10) 

The  functions  Ji{U)  is  a  Bessd  functions  of  argument  U  and  order  1  and  Hi{U)  is 
a  Struve  function  of  argument  U  and  order  1  where: 

£r  =  *on,-(i'  +  i).  (11) 


For  dielectric  meteriali,  the  diffraction  coefficient  including  caustic  correction  is  pven 
by: 


where  7«(£^)  is  as  described  previously.  Equation  (12)  works  well  in  describing  the 
far  sone  fields  from  a  disk  where  a  flat  face  is  always  visible  and  will  later  be  used 
extensively  in  analyzing  the  scattering  from  trees. 


SECTION  3 
Background 

1  Introduction 

Before  the  acattenng  from  « tree  can  be  analysed,  some  background  information  on 
tree  modeling  is  needed.  In  this  Chapter,  a  summary  of  current  methods  of  modeling 
a  tree  is  presented.  The  summary  includes  a  description  of  what  tree  components 
need  to  be  consideted  and  how  these  components  are  usually  modeled.  Finally, 
information  is  provided  on  modeling  the  ground. 

2  Modeling  of  a  lYee 

Many  different  methods  are  used  to  model  a  tree  depending  on  the  species  of  the 
tree  to  be  modeled  and  the  frequency  range  of  interest.  For  the  purposes  of  this 
study,  an  attempt  will  be  made  to  use  average  values  for  the  tree  parameters  such 
that  a  wide  range  of  tree  types  are  represented.  The  primary  frequency  range  of 
interest  is  between  0.3-L3GHz  although  other  frequencies  are  studied. 

Due  to  the  low  frequencies  of  interest,  amplifications  can  be  made  in  modeling  the 
trunk  and  branches  of  the  tree.  For  the  frequencies  between  0.3*1.3GHa,  the  surface 
of  the  tree  and  branches  can  be  treated  as  smooth  rince  the  rou^mess  of  the  tree  bark 
win  be  smaU  compared  to  a  wavelength.  Also,  the  leaves  and  smafier  branches  on 
the  tree  are  ignored  when  modeling  the  trunk  and  branches.  However,  for  a  complete 
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modd  of  a  tree,  the  leaves  and  small  branches  need  to  be  accounted  for.  The  only 
remaining  features  of  the  tree  are  the  trunk  and  the  larger  brandies.  The  trunk  and 
brandies  axe  then  modeled  as  dielectric  cylinders  with  a  single  finite  length  cylinder 
representing  the  trunk  and  smaller  finite  length  cylinders  representing  the  larger 
branches.  There  axe  many  dififerent  methods  used  to  pick  the  height,  diameters, 
and  dielectric  properties  of  the  cylinders.  The  different  techniques  of  picking  these 
parameters  depend  laxgdy  on  the  tree  spedes  and  what  seems  to  work. 

The  choice  of  cylinder  sues  used  to  modd  the  tree  trunk  has  been  reported  by 
maxiy  authors.  Specifically,  formulas  by  M.  0.  Kolawole  [12],  and  G.  Sun  and  D. 
S.  Simonett  [13]  axe  summarised  here.  Ki^wole  uses  a  simple  relationship  between 
the  diameter  of  the  tree  at  breast  height  {DBH)  to  the  tree  height  h  and  effective 
he%ht  kt  pven  by: 

DBH  =  0.00532  +  0.0172h  (13) 

(14) 

where  all  dimensions  axe  in  meters.  G.  Sun  and  D.  S.  Simonett  use  another  rela^ 
tkmship  between  the  DBH  of  a  tree  and  the  tree  height  given  by: 

DBH  =  0.00532  +  0.0im  (15) 

K  -  0.59ih-.82  (16) 

ht  =  hv  +  -{h  —  hy)  (17) 

where  again  all  dimensions  axe  in  meters.  In  Equation  (17),  h«  is  the  distance  from 
the  ground  to  the  widest  part  of  the  crown.  The  parameter  ht  is  the  length  of  the 
dielectric  cj^der  that  is  used  to  modd  the  trunk  of  the  tree,  i.e.,  the  effective 
hd|^t  of  the  tree.  Usually,  the  diameter  of  the  tree  {DBH)  is  known  and  is  used  to 
solve  for  the  heights  and  hi.  These  relationships  between  the  tree  diameter  and 
effective  height  axe  used  repeatedly  in  modding  tree  trunks. 

While  there  is  infonnation  on  the  dimenuons  of  the  tree  trunk  used  for  modeling, 
there  is  little  discusnon  of  the  dimensions  of  branches.  The  fiteratnxe  is  short  on 
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detailf  for  the  lengths  and  diameters  that  are  used  to  model  tree  branches.  Therefore, 
an  intuitnre  approach  for  selecting  diameters  and  lengths  will  be  used  based  upon  the 
exact  modeling  needs  when  analysing  their  scattering.  As  a  result,  the  discussion 
of  the  lengths  and  diameters  of  branches  for  modeling  will  be  postponed  to  the 
discussion  of  branch  scattering. 

3  Dielectric  Properties  of  T^ees 

In  addition  to  knowing  the  siae  of  a  tree,  its  dielectric  properties  must  also  be  known 
for  realistic  modeling.  There  have  been  many  papers  that  present  measured  dielec¬ 
tric  constants  of  trees  [14,  15].  The  daily  and  seasonal  variations  of  the  dielectric 
constants  of  trees  have  also  been  examined  by  several  authors  and  for  different  tree 
species  [14, 16].  But,  for  the  purposes  of  this  study,  the  variation  in  tree  parameters 
win  be  ignored  and  instead  a  set  of  values  for  the  dielectric  constants  wiU  be  used 
that  should  cover  ty]ucal  values  including  the  variations  in  the  dielectric  constant. 
In  theory,  however,  the  variations  can  be  taken  into  account. 

The  modeling  of  the  dielectric  properties  of  wood  can  be  divided  into  two  groups, 
the  use  of  a  single  homogeneous  piece  of  dielectric  and  the  use  of  multiple  layers  of 
dielectric  to  rimulate  the  individual  parts  of  the  tree  structure.  While  the  simplicity 
of  modeling  the  wood  as  a  tingle  la;^  of  dielectric  is  apparent,  it  can  be  argued 
that  the  multilayer  approach  is  more  accurate  since  it  is  more  faithful  to  the  actual 
structure.  This  is  because  the  dielectric  constant  is  lower  on  the  inside  of  a  tree. 
The  dielectric  constant  then  increase  as  a  function  of  radius  before  decreasing  to  a 
lower  ditiectric  constant  near  the  bark.  This  profile  of  the  dielectric  constant  forms 
the  basis  of  the  3-layer  modding.  Both  methods  of  modeling  a  tree  as  tin|^  and 
multiide  layers  are  presented  here.  In  Table  1,  the  relative  didectric  constants  of 
wood  used  in  tingle  layec  modeling  are  presented.  The  wood  is  assumed  to  have  no 
magnetic  properties  such  that  the  permeability  of  the  wood  pe  is  equal  to  /iq.  The 
wide  range  of  values  in  l^ble  1  shows  the  large  differences  in  modding  trees  with 
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IkUe  1:  Dielectric  constants  used  to  model  the  dielectric  properties  of  wood  at 
L-Band  from  the  literature. 


«r 

author(s)  and  reference 

5-jl.5 

F.  Amir,  et.  al.  [17] 

12.4-j4.9 

M.A.  Karam  and  A.K.  Fung  [18] 

18-j6 

5-jl.7 

24-j8 

Kyle  McDonald  [19] 

4-jl 

Tsan  Mo  and  J.R.  Wang  [20] 

20-j5 

Guo-qing  Sun,  et.  al.  [21] 

Table  2:  Didectric  constants  used  to  model  the  dielectric  properties  of  wood  at 
L-Band  from  the  literature  for  multilayered  modeling. 


layer 

layer  thickness  (cm) 

m 

inner  layer 

10 

15-j7 

outer  layer 

0.5 

4-jl 

dielectric  materials.  The  large  'variation  in  dielectric  constants  is  explained  by  the 
wide  range  of  tree  spedes  that  have  been  measured.  The  variation  in  Table  1  is  also 
due  to  the  seasonal  and  daily  variations  in  the  dielectric  constants  of  wood  mainly 
due  to  the  diange  in  moisture  content  of  the  wood.  In  Table  2,  the  rdative  dielectric 
e<mstants  and  layer  information  are  presented  for  2-layer  modeling.  In  addition  to 
using  a  2-layer  model,  a  third  layer  can  be  added  to  to  simulate  more  accurately  the 
didectric  profile  of  the  tree  as  mentioned  previously.  The  use  of  3-layer  modding  is 
discussed  in  later  chapters.  The  differences  in  modeling  the  cylinders  as  singile  and 
multiple  layers  are  examined  further  in  following  chapters. 
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4  Dielectric  Properties  of  Soil 

The  effect  of  the  ground  on  the  scattering  from  a  tree  can  be  of  importance.  There¬ 
fore,  a  model  for  the  ground  needs  to  be  developed.  Typically,  for  the  primary 
frequencies  of  interest,  the  wavelength  u  long  enough  so  that  the  ground  can  be 
treated  as  smooth.  The  ground  is  also  assumed  planar  near  the  tree.  The  sc^  is 
then  usually  modeled  as  a  homogeneous  dielectric  half-space.  Table  3  summarises 
typical  dielectric  properties  for  sdl  that  will  be  used  in  modeling  the  scattering  from 
a  tree. 

Table  3:  Dielectric  constants  used  to  model  the  dielectric  properties  of  soil  at  L-Band 
from  the  literature. 


<r 

comments  and  source 

7-jl.5 

moist  soil  foom  hiichael  0.  Kolawole  [12] 

18-j5 

wet  sdl  from  Michael  0.  Kolawole  [12] 

20-j0.15 

soil  at  Ohio  State  University  Radio  Telescope  [22] 

SECTION  4 

Scattering  From  Logs 

1  Introduction 

A  model  for  the  icatteruig  characteristics  of  wood  needs  to  be  developed  that  can 
be  used  in  the  analysis  of  the  scattering  from  trunks  and  branches.  In  Chapter  3,  a 
summary  d  the  parameters  used  in  modeling  the  various  parts  of  a  tree  is  presented 
along  with  dielectric  constants.  However,  due  to  the  wide  range  of  values  that  are 
presented,  it  is  desired  to  find  which  values  and  techniques  are  appropriate  for  the 
cases  of  interest  in  this,  analysis.  To  do  this,  the  parameters  and  ideas  that  are 
presented  in  Chapter  3  are  used  to  develop  a  model  that  is  appropriate  to  the  cases 
of  interest  in  this  analysis. 

To  develop  this  model,  the  scattering  from  a  log  is  measured  and  theory  is  used 
to  match  the  measurements.  In  doing  this,  a  set  of  parameters  fi>r  the  wood  is 
found  that  will  allow  theory  to  predict  the  scattering  frmn  the  log.  These  wood 
parameters  will  then  be  used  later  to  model  the  wood  in  more  comidex  geometries 
invdvittg  branches  and  trunks. 

To  find  the  parameters  for  the  wood,  a  method  of  modding  the  scattering  from 
a  log  is  presented.  First,  methods  of  accounting  for  the  effects  of  log  roughness  and 
dielectric  constant  on  the  calculated  scattered  fidd  of  the  log  are  discussed.  Then 
measured  results  for  the  RCS  of  three  logs  are  presented.  The  methods  of  modeling 
a  log  using  the  roughness  and  dielectric  constants  are  then  applied  to  compare  to 
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the  measured  RCS  of  the  logs.  The  resulting  values  for  the  dielectric  constants  and 
roughness  used  to  match  the  measured  log  RCS  provides  a  set  of  values  for  the  wood 
that  can  be  used  to  model  the  characteristics  of  wood  in  branches  and  trunks. 

2  Theory  of  Modeling  a  Log 

To  model  the  scattering  from  a  log,  three  things  are  taken  into  account.  First, 
a  modd  for  the  dielectric  properties  of  the  wood  must  be  found  that  accurately 
represents  the  dielectric  constants  of  the  wood  as  a  function  of  frequency.  Then,  a 
method  of  accounting  for  the  roughness  of  the  log’s  bark  must  be  chosen  depending 
on  the  frequencies  of  interest.  Finally,  the  dielectric  and  roughness  information  must 
be  included  into  a  solution  for  the  scattering  for  a  finite  length,  rough,  dielectric 
cylinder. 

2.1  Dielectric  Properties  of  Wood 

To  model  the  dielectric  properties  of  wood,  the  variation  of  the  dielectric  constant 
with  frequency  must  be  taken  into  account.  Various  didectric  constants  have  been 
published  to  modd  wood  and  are  summarized  in  Table  1.  These  values  represent 
tyjncal  values  used  in  L-band.  However,  a  method  of  calculating  the  scattered  fidd 
over  a  broad  range  of  frequendes  u  of  interest.  There  are  many  papers  that  present 
data  on  the  dielectric  constant  of  wood  in  specific  frequent  bands  but  there  is 
little  data  available  that  present  the  dielectric  constant  as  a  function  of  frequency. 
Therefore,  a  method  will  be  used  to  vary  the  didectric  constant  as  a  function  of 
frequency  using  the  values  as  presented  in  Table  1. 

Bar  the  purposes  of  this  analyns,  the  real  part  of  the  didectric  constant  will 
always  be  hdd  constant  with  frequency.  This  is  done  to  simplify  the  analysis  and 
is  oonastent  with  the  small  variation  in  the  real  part  of  the  didectric  constant  with 
frequency  shown  in  the  literature. 


Although  the  real  part  of  the  dielectric  constant  varies  little,  the  imaginary  part 
has  been  shown  to  vary  considerably  with  &equency[14].  Therefore,  a  simple  method 
of  varying  the  ima|^ary  part  of  the  dielectric  constant  with  frequency  will  be  dis¬ 
cussed  and  compared  with  holding  the  dielectric  constant  fixed  with  frequency.  Tak¬ 
ing  the  values  from  Table  1  and  by  knowing  the  approximate  frequency  at  which 
the  dielectric  constants  were  measured,  the  conductivity  of  the  material  at  that 
frequency  can  be  found  by: 

fft  =  •^WMip[£r]  X  UtfO  (18) 

where  Ut  is  the  frequency  at  which  the  relative  dielectric  constant  Cr  is  measured. 
The  value  of  a,  from  Equation  (18)  is  then  held  constant  as  the  frequency  is  varied. 
The  relative  dielectric  constant  6r  then  varies  with  frequency  as: 

€r  =  — (e  -h  ^)  wh«c  tTe  “id  e  are  fixed.  (19) 

Co  JIt) 

This  method  of  varying  the  dielectric  constant  will  be  corapared  to  the  method  of 
holding  the  dielectric  constant  fixed  with  frequency  later. 

2.2  Effect  of  Wood  Roughness 

The  roughness  of  the  log  bark  also  must  be  taken  into  account.  While  the  roughness 
of  the  bark  will  be  small,  as  the  frequenqr  increases  the  effect  of  the  roughness  will 
increase.  For  the  logs  and  frequendes  of  interest,  an  appropriate  model  for  treating 
the  roughness  is  used  where  the  roughness  factor  is  a  function  of  the  surface  hdght 
variability  [23].  This  model  is  only  valid  for  slightly  rough  surfaces  where  the  main 
reflection  from  the  surface  is  the  specular  component  and  the  diffuse  scatter  caused 
by  the  surface  roughness  is  small.  The  modd  modifies  the  r^ection  coefildent  by: 

Trough  =  ^.r  (20) 

where: 

{SrY  =  (21) 
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(22) 


AirAk  . 
—r—sm'i. 


In  Equations  (20)-  (22),  F  is  the  reflection  coefficient  from  the  smooth  surface,  Trough 
is  the  modified  reflection  coeffident  for  the  rough  surface.  Ah  is  the  maximum  hdght 
difference  between  two  points  on  the  rough  surface,  and  7  is  the  angle  between  the 
inddent  wave  and  the  surface.  This  formula  is  valid  only  when  the  height  variation 
Ak  is  less  than  which  is  true  for  the  cases  of  interest.  Values  of  Ah  will  depend 
on  the  type  of  wood  and  number  of  bumps  in  the  wood  bark  and  will  be  found  later 
when  the  wood  measurements  are  modeled. 


2.3  Solution  of  Dielectric  Cylinder  Scattering 

The  scattering  from  a  finite  length  dielectric  cylinder  must  be  calculated.  To  ac¬ 
complish  this,  the  solution  for  scattering  by  a  multilayered  infinite  length  dielectric 
cylinder  described  in  Appendix  B  is  used.  Since  this  solution  is  for  infinite  length 
cylinders  it  must  be  modified  for  finite  length  cylinders.  A  simple  method  of  doing 
this  is  to  convert  the  2-dimensional  echo  width  to  the  S-dimensional  radar  cross 
section  (RCS).  This  is  accomplished  by  the  well  known  formula: 

2X* 

osD  *  —  X  a2D  (23) 

where  L  is  the  length  of  the  finite  length  cylinder.  This  conversion  pves  accurate 
answers  when  the  length  L  is  much  greater  than  the  cylinder  diameter  and  A. 


3  Measurement  of  the  Radar  Cross  Section  of 
Logs 

To  measure  the  radar  cross  section  of  a  log,  a  tree  apprmdmately  4”  (0.10m)  in 
diameter  was  cut  down  from  the  woods  near  the  Ohio  State  University  Big  Ear 
radio  tdescope  where  forest  propagation  measurements  have  been  made  in  the  past. 
The  tree  was  a  sugar  maple  which  comprise  a  large  number  of  the  trees  in  the 
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woods.  The  tree  trunk  was  cut  into  3  lengths  approximately  24”  (0.61m)  long.  The 
diameters  of  the  logs  varied  from  approximately  3.5”-4”  (0.089-0.10m).  The  bark 
of  the  logs  was  fairly  smooth  with  main  groves  in  the  bark  approximately  2-3mm 
deep.  The  branches  were  cut  off  leaving  a  few  small  bumps  on  each  log.  After 
the  tree  was  cut  into  logs,  the  ends  and  the  bark  were  spray  coated  with  shellac  to 
reduce  moisture  loss  due  to  evaporation.  Although  the  layer  of  shellac  has  a  slightly 
different  dielectric  constant  than  the  wood,  the  layer  of  shellac  was  very  thin  and, 
therefore,  is  thought  to  have  negfipble  impact  on  the  wood  dielectric  properties. 
Also,  after  the  Ic^  were  removed  from  the  woods,  they  were  wrapped  in  plastic 
bi^s  to  further  reduce  water  loss. 

The  logs  were  then  measured  in  a  compact  range  the  next  day  &om  2-18GHz.  The 
logs  were  measured  in  this  frequency  range  due  to  necessity  of  available  equipment 
at  the  time  of  the  measurements.  Although  these  frequencies  are  higher  than  the 
primary  frequencies  of  interest,  the  results  obtained  by  the  measurements  should 
be  applicable  to  the  frequencies  between  0.3>1.3GHz.  The  measurements  consisted 
of  specular  impulse  responses  from  the  smoothest  section  of  the  logs  (i.e.,  the  side 
with  the  fewest  bumps).  The  measurements  were  then  repeated  19  days  later  to 
allow  for  some  drying  of  the  logs  to  see  how  a  small  decrease  in  the  water  content  of 
the  logs  would  affect  the  RCS.  For  the  second  set  of  measurements,  the  orientation 
of  the  log  was  maintained  as  dose  as  possible  to  the  ori^al  measurements.  Both 
dayu  of  measurements  are  combined  and  are  shown  for  the  three  logs  in  Figures  2, 3, 
and  4.  The  RCS  is  represented  in  dedbels  with  respect  to  a  square  meter  for  the 
transverse  magnetic  (TM)  case  with  the  electric  field  aligned  with  the  log  mtik  and 
the  transverse  dectric  (TE)  case  with  the  magnetic  fidd  aligned  with  the  log  axis. 

It  is  seen  from  the  measurements,  that  the  RCS  of  each  log  didn’t  change  ap- 
predably  fior  the  2  days  of  measurements  at  the  frequendes  bdow  approximatdy 
lOGHz.  But  above  lOGHz,  the  RCS  for  the  2  days  of  measurements  for  each  log 
tended  to  dian^.  This  variation  can  be  attributed  to  the  difficulty  in  predsdy  align- 
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Figure  2:  Measured  RCS  of  log  #1  for  TM  aad  TE  pcdaiiaations. 
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f^guze  3:  Measiized  RCS  of  log  #2  for  TM  and  T£  penalizations. 
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Fignte  4:  Measiixed  RCS  of  U>g  #3  for  TM  and  TE  polaxisations. 
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isg  the  log  to  the  same  location  as  the  pxevions  days  measurement  since  a  rotation 
abont  the  long  axu  of  the  log  was  difficult  to  contrcd.  Since  the  logs  don’t  have  a 
perfectly  circular  cross  and  are  not  entirely  straight,  a  small  rotation  causes  a  change 
in  the  surface  characteristics  of  the  log  being  measured.  Also,  the  resulting  small 
rotation  of  the  log  caused  a  different  part  of  the  log  surface  to  be  measured,  thereby, 
the  effective  surface  roughness  changed  slightly.  Since  the  surface  roughness  of  the 
log  has  an  larger  effect  on  RCS  at  higher  frequencies  as  shown  by  Equation  (22),  the 
small  rotational  error  caused  a  large  difference  on  the  RCS  at  higher  frequencies.  By 
the  same  reasoning,  the  small  differences  in  the  data  for  the  2  days  of  measurements 
at  the  lower  frequencies  show  that  the  surface  roughness  has  a  small  effect  on  the 
RCS  of  the  logs  at  the  lower  frequencies  and  that  any  drying  of  the  logs  also  had  a 
small  effect. 


4  Comparison  of  Log  RCS  with  Theory 

It  is  now  desirable  to  try  to  relate  the  different  techniques  of  modeling  logs  from 
Section  2  to  the  measurements  of  Section  3.  To  accomplish  this,  the  surface  rough- 
new  model  described  in  Equations  (22)  and  (20)  is  used  to  modify  the  echo  width  in 
the  exact  solution  pven  by  Equation  (76).  The  scattered  field  in  Equation  (76)  was 
multipfied  by  the  coeffident  S,,  This  results  in  the  echo  width  for  a  rou^  cylinder 
^tD  rcttgh  where; 

2A  **  * 

V2D-rm>^-  -  S,'^enDnCOtn^  .  (24) 

nsO 

The  edio  width  given  by  Equation  (24)  is  then  modified  from  the  two-dimensional 
(2D)  echo  width  to  the  three-dimensional  (3D)  RCS  by  Equation  (23).  Dielectric 
constants  are  then  chosen  from  Thble  1  and  are  used  to  model  the  log  as  a  sinf^e- 
layer  homogeneous  dielectric  cylinder.  The  length  of  log  #2  is  27.5”  (0.699m)  and 
an  average  value  for  the  diameter  for  log  #2  of  3|"  (0.092m)  is  chosen  and  used  to 
model  the  log.  Two  cases  are  run,  one  where  the  conductivity  is  held  fixed  and 
the  other  where  the  dielectric  constant  e.  is  held  fixed.  The  two  cases  are  shown  in 


Figure  5  for  er»24-j8,  and  for  Rieal[er]=24  and  «re=0.5563.  Figure  5  iliows  that  the 
two  curves  follow  the  same  general  trend,  however,  the  constant  curve  shows  a 
stronger  creeping  wave  component  than  the  constant  Cr  curve.  The  creeping  wave 
component  can  travel  ea^y  around  the  smooth  cylinder  used  in  the  solution.  But, 
since  the  roughness  of  tree  bark  will  quickly  attenuate  the  creeping  wave  component, 
the  creeping  wave  component  of  the  sdution  must  be  removed  by  processing  the 
scattered  field.  The  creeping  wave  is  removed  by  processing  the  impulse  response  of 
the  theoretical  solution  and  filtering  out  the  creeping  wave.  The  processed  results 
for  the  same  cylinders  as  calculated  in  Figure  5  are  shown  in  Figure  6.  By  examining 
Figure  6,  once  the  creeping  wave  component  is  removed  from  the  scattered  fidd,  it 
doesnH  matter  whether  <r«  or  Cr  is  held  constant  with  frequency  since  the  results  are 
the  same  once  the  cree|»ng  wave  is  removed.  Therefore,  since  both  methods  produce 
the  same  result  after  the  creeping  wave  component  is  removed,  future  calculations 
will  keep  the  dieJectric  constant  e,  constant  with  frequent  when  the  creeping  wave 
component  is  removed. 

Since  the  three  logs  came  from  the  same  tree  trunk,  the  measured  RCS  can  help 
illustrate  the  variations  in  the  RCS  of  a  tree.  To  show  this  variation,  the  RCS  of 
the  three  logs  is  converted  to  the  2D  edio  width  by  equation  (23).  This  is  done 
to  remove  the  effect  of  the  different  log  lengths.  The  resulting  echo  widths  for  the 
three  logs  and  2  different  days  of  measurements  is  shown  in  Figures  7  and  8  for  TE 
and  TM  polarisations,  respective^.  The  purpose  of  these  plots  is  not  to  show  the 
individual  edbo  width  for  each  log  measurement  but  to  show  the  variation  in  the 
scattermg  from  the  logs.  At  the  low  frequencies  near  2.0GHs,  the  echo  width  varies 
about  4dB  while  at  higher  frequencies  the  fluctuation  is  almost  15dB.  This  shows 
that  for  the  same  tree  trunk,  the  RCS  of  different  sections  will  vary  by  several  dB. 
Due  to  this  fluctuation,  it  is  difficult  to  model  the  RCS  of  a  log  (or  tree  trunk)  with 
a  high  degree  of  accuracy.  Therefore,  it  can  be  «q>ected  that  any  modeling  of  a  tree 
trunk  will  only  be  accurate  to  within  a  few  dB. 
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FSgaxe  5:  Compaxuon  of  the  RCS  of  the  dielectric  cylinder  used  to  model  a  log  for 
relative  dielectric  cmutaati  ers24>j8  and  Realfcr]— 24,  0’e=O.5563. 
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P^^nze  7:  Echo  width  for  TE  p<daxisation  of  logs  #1,  #2,  and  #3  from  processing 
RCS  measurements  to  show  variation  in  echo  width  for  logs  from  the  same  tree. 
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Echo  Width  (dBm) 


Frequency  (GHz) 

Flgue  8:  Echo  wwidth  for  TM  polarigation  of  logs  #1,  #2,  and  #3  from  processing 
RCS  measurements  to  show  wariatirm  in  echo  width  for  logs  from  the  same  tree. 
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The  rooghneti  modifiattion  for  clighily  rough  surfaces  as  described  in  Equa- 
timis  (22)  is  then  added  to  the  modeling  of  the  logs.  In  Figures  9  and  10,  the  echo 
width  for  the  3  log  measurements  are  compared  with  two  different  dielectric  models 
for  the  1<^.  The  models  are  for  dielectric  constants  of  er=12.4-j4.9  and  er=24-j8 
with  toughness  factor  (Ah)  of  2mm  and  3mm  since  the  log  roughness  is  between 
2mm  and  3mm.  The  models  for  the  logs  are  within  approximately  2dB  lower  than 
the  measured  echo  width  at  the  lower  frequencies  and  are  within  the  range  of  mea¬ 
surements  at  the  higher  frequencies.  Since  it  is  difficult  to  modd  the  logs  more 
accurately,  these  dielectric  constants  are  considered  adequate  to  represent  the  range 
of  dielectric  properties  for  the  trees  of  interest  in  modeling.  Based  on  these  results, 
for  future  calculations  of  the  RCS  from  wood,  the  dielectric  constants  used  to  model 
the  wood  win  either  be  er=12.4-j4.9  or  Cr=s24-j8.  Since  using  roughness  factors  (Ah) 
of  2-3mm  falls  within  the  measured  results  at  higher  frequencies,  when  the  roughness 
is  important  the  value  of  Ah  wiQ  be  between  2mm  and  3mm.  These  values  should 
provide  accurate  modeb  for  actual  wood. 

The  sur&ce  roughness  has  a  significant  effect  on  the  RCS  of  the  logs  at  higher 
frequencies  while  at  lower  frequencies,  the  surface  roughness  has  a  smaller  effect.  In 
Figures  9  and  10,  a  change  in  Ah  from  2-3mm  causes  a  large  change  in  the  echo  width 
of  about  lOdB  at  18GHz.  This  sensitivity  u  farther  iUustrated  in  Figure  11  which 
shows  the  roughness  effect  for  various  values  of  Ah.  It  is  seen  that  the  roughness  only 
has  an  effect  of  about  3dB  at  6GHs  for  a  Ah=3mm.  Therefore,  when  the  roughness 
factor  is  removed  from  the  echo  width  calculations,  the  difference  between  the  theory 
and  the  measurements  will  be  within  the  variation  of  the  measured  echo  width  for  the 
logs  of  4dB  at  6GHs.  Therefore,  the  roughness  frctor  can  be  ignored  at  frequencies 
bdbw  6GHs  and  the  log  surface  can  be  treated  as  smooth  without  causing  significant 
error.  For  the  primary  frequencies  of  interest  between  0.3-1.3GHz,  the  effect  of  the 
surface  roughness  is  even  smaller  and  therefore,  it  can  be  ignored  as  was  discussed 
in  Chapter  3.  Above  6GHx,  the  sensitivity  to  the  roughness  factor  is  too  great  to  be 
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Hguxe  9:  Echo  width  for  TE  poUxinition  of  logs  #1,  #2,  and  #3  from  processing 
RCS  measurements  compared  to  different  didectxic  models  for  the  logs. 
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Figuie  10:  Echo  width  for  TM  polarization  of  logs  #1,  #2,  and  #3  from  processing 
RCS  measurements  compared  to  different  dielectric  models  for  the  logs. 
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used  in  modeling  and,  therefore,  is  not  recommended  since  it  is  difficult  to  measure 
the  Ah  value  and  then  calculate  the  RCS.  Instead  the  value  of  Ah  must  be  chosen 
to  match  the  measurements.  Therefore,  this  technique  of  modding  the  roughness  of 
the  log  is  not  recommended  at  high  frequencies. 


Figure  11:  Variation  of  roughness  factor  5,  for  different  values  of  Ah  as  a  function 
of  frequency. 

However,  if  the  roughness  factor  is  varied  until  the  calculated  RCS  matdi  the 
measured  RCS,  the  modeling  of  the  log  closely  approximates  the  measured  RCS. 
In  Figure  12,  comparison  between  the  measured  and  calculated  RCS  for  log  #2  is 
shown.  For  these  plots,  the  dielectric  constant  is  dbosen  as  6r=24-j8.  The  surface 
lou^mess  is  found  by  a^insting  Ah  until  the  calculated  RCS  matched  the  measured 
data.  While  this  method  seems  rather  unsdentiiic,  for  the  log  geometry,  there  is  no 
way  to  actually  calculate  the  roughness  factor  Ah  directly  due  to  the  complexity  of 
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Figure  12:  Comparison  between  the  measured  and  calculated  RCS  of  log  #2.  For 
the  modeling,  the  parameters  are  er=24-j8  and  Ah=0.0025m 
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the  log  fuzface.  Therefore,  once  a  value  is  found  by  trial  and  error,  it  is  verified  that 
this  value  is  reasonable  for  the  log.  A  roughness  Ah=0.0025m  is  needed  to  match 
the  measured  data.  This  roughness  is  on  the  same  level  of  roughness  that  can  be 
measured  for  the  bark  of  between  2*3mm  which  doesn’t  include  large  features  such 
as  bumps  from  branches.  Therefore,  it  is  felt  that  this  is  a  reasonable  value  for  the 
roughness  factor  Ah.  Also,  for  Figure  12  the  creeping  wave  component  is  processed 
out  as  discussed  previously.  So,  this  method  of  accounting  for  the  surface  roughness 
can  provide  good  agreement  with  measurements.  However,  as  stated  earlier,  since 
the  RCS  from  modeling  the  log  is  very  dependent  on  the  surface  roughness  Ah,  and 
there  is  no  way  to  measure  Ah  directly,  this  method  of  modeling  the  logs  at  higher 
frequencies  should  be  used  with  caution. 

5  Summary 

Methods  of  calculating  the  effect  of  the  dielectric  constants  and  roughness  of  a  log 
on  the  calculated  scattered  field  is  presented.  The  measured  RCS  of  three  logs 
is  then  shown.  By  comparing  the  theory  of  calculating  the  RCS  of  a  log  to  the 
measurements,  the  effect  of  roughness  and  didectric  constant  is  observed.  It  is 
found  that  the  roughness  of  the  bark  had  a  significant  contribution  to  the  RCS 
of  the  logs  above  approodmatdy  6GHz  but  had  a  much  smaller  impact  at  lower 
frequendes.  Therefore,  for  future  modeling,  the  surface  roughness  can  be  ignored  for 
frequendes  bdow  6GHz.  The  effect  on  varying  the  didectric  constant  with  frequency 
versus  holding  it  fixed  is  also  examined.  It  is  found  that  when  holding  the  didectric 
constant  fixed,  the  creeping  wave  component  decreased  with  frequency  as  compared 
to  the  case  where  the  dielectric  constant  varies  with  frequency  in  which  case  the 
creeping  wave  remained  significant  at  hi^er  frequendes.  However,  the  method  used 
to  calculate  the  scattered  fidd  was  for  a  smooth  cylinder  and,  therefore,  a  creeping 
wave  could  be  supported.  But  for  an  accurate  modeling  of  the  log  to  fully  account 
for  the  surface  roughness,  the  creejung  wave  needs  to  be  removed  since  the  rough  log 
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fuxfmce  reduces  tlie  creeping  waves.  When  the  creeping  wave  is  removed,  the  results 
for  a  constant  and  a  varying  dielectric  constant  are  identical.  TherdTore,  for  future 
modeling,  the  dielectric  constant  can  be  chosen  to  vary  or  remain  fixed  and  will 
produce  the  same  results  as  long  as  there  is  no  creeping  wave  components.  Finally, 
the  effect  of  choice  of  the  dielectric  constants  is  examined.  It  is  found  that  using 
dielectric  constants  of  er=24-j8  or  £r=12.4-j4.9  for  the  wood  produced  the  correct 
scattered  field  magnitude  within  approximately  2dB  of  the  measured  values.  The 
difference  of  2dB  is  shown  to  be  small  and  a  good  enough  match  since  the  echo 
widths  of  the  three  logs  vary  by  more  than  4dB  for  logs  from  the  same  tree  trunk. 
Therefore,  these  values  for  the  wood  dielectric  constants  will  be  used  extensively  to 
model  wood  in  later  chapters. 
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SECTION  5 

Scattered  Fields  in  the  Vicinity  of 
a  TVee  TVunk 

1  Introduction 

The  scattered  fields  in  the  vicinity  of  a  cylinder  used  to  model  a  tree  trunk  are  of 
interest.  In  Chapter  4,  the  modeling  of  a  log  as  a  dielectric  cylinder  is  discussed  and 
agreement  between  the  modeling  and  measurements  is  shown.  This  chapter  discusses 
the  scattered  field  in  the  vicmity  of  the  dielectric  cylinder  used  to  represent  the  tree 
trunk.  The  dielectric  parameters  used  for  the  dielectric  cylinder  are  based  in  part  on 
the  results  of  the  previous  chapter  and  also  on  the  modding  discussed  in  Section  2. 
The  purpose  of  the  analysis  is  tc^  look  for  trends  in  the  near-fidd  scattering  of  a 
cylinder  for  various  geometries. 

The  analysis  is  comprised  of  two  parts.  The  first  part  involves  studying  the 
near  fidd  scattering  of  a  cylinder  representing  a  tree  trunk  using  an  exact  numeri¬ 
cal  solution  of  a  multilayered  didectric  cylinder  of  infinite  length  under  plane  wave 
inddence.  Various  combinations  of  qrlinder  parameters  are  used  to  discover  trends 
in  the  scattered  patterns  for  various  cylinder  geometries.  The  second  part  at  the 
aaaljrsis  invdves  comparing  the  applicabifity  of  the  above  exact  solution  for  infinite 
length  cylinders  with  normal  inddence  to  finite  length  cylinders  with  oblique  ind¬ 
dence.  This  is  done  to  find  a  range  of  angles  over  which  the  solution  for  an  infinite 
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lengtli  cjlinder  with  normal  incidence  can  be  used  for  a  finite  length  cylinder  with 
obfiqne  incidence.  For  this  comparison,  a  UTD  solution  for  a  finite  length  perfect 
electrical  conducting  (PEC)  cylinder  alone  and  above  an  infinite  ground  plane  with 
oblique  plane  wave  incidence  is  compared  to  the  above  exact  solution  for  an  infinite 
length  PEC  cylinder.  A  UTD  solution  for  PEC  cylinders  is  used  since  a  more  gen¬ 
eral  solution  for  the  near  field  of  a  multilayered  dielectric  cylinders  of  finite  length 
and  oblique  incidence  is  not  presently  available.  From  this  comparison,  a  range  of 
applicability  is  found  which  should  be  valid  for  the  multilayered  didectric  case. 

2  Target  Geometries 

The  geometry  used  for  the  analysis  of  infinite  length  cylinders  is  shown  in  Figure  13. 
It  consists  of  a  multilayered  dielectric  blinder  of  infinite  length  with  circular  cross 
section  in  the  z-y  plane.  A  plane  wave  is  incident  from  the  n^ative  z-axis.  The 
solution  to  this  geometry  is  discussed  in  Appendix  B. 

The  geometry  used  in  the  comparison  cases  for  oblique  incidence  is  shown  in  Fig¬ 
ure  14.  The  UTD  solution  for  a  finite  length  PEC  cylinder  with  oblique  plane  wave 
incidence  is  calculated  using  the  Numerical  Electromagnetic  Code  -  Basic  Scattering 
Code  (NEC-BSC)  code  [24].  This  code  does  not  accuiatdy  provide  the  scattered 
fidds  in  the  deep  shadow  region  of  a  didectric  cylinder.  Therefore,  only  the  PEC 
case  is  analysed  for  the  finite  length  qrlinder.  The  NEC-BSC  code  outputs  the  total 
fidd  for  a  radial  180  d^ree  cut  in  the  forward  scatter  repon  at  a  constant  radial 
distance  R  about  the  center  of  the  finite  length  cylinder  with  the  angle  ^  defined  as 
in  figure  13. 

The  geometry  fm  the  case  of  the  finite  length  cylinder  above  an  infinite  ground 
plane  used  to  simulate  the  ground  is  shown  in  Figure  15.  This  geometry  is  analjrzed 
since  it*s  doser  to  the  actual  geometry  that  would  occur  in  nature.  This  case  pro¬ 
duces  the  same  type  of  output  as  the  case  without  the  ground  plane  but  allows  the 
height  dl  the  recdver  above  the  ground  plane  to  be  varied. 
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Figure  14:  Geometry  used  for  analytis  of  fimte  length  cylinders  with  oblique  plane 
wvre  inddenoe. 


Figure  15:  Geometry  used  for  analyns  of  finite  length  cylinders  with  oblique  plane 
wfcve  madenoe  inddence  above  an  infinite  ground  plane  to  simulate  the  ground. 
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3  Analysis  and  Results 

For  the  analyab,  different  combinations  of  material  properties  and  thicknesses  for  the 
cylinder  are  used  to  find  any  trends  in  the  near-sone  scattered  field  to  cylinder  model 
parameters.  The  dielectric  constants  ^  and  dimensions  of  the  dielectric  cylinder  nsed 
are  ficom  Tables  1  and  2.  Chapter  4  shows  that  the  single  layer  dielectric  modd  is 
snfiident  to  match  with  measurements  for  the  backscatter  RCS.  However,  nnce  the 
fields  in  the  vidnity  of  the  cylinder  are  now  of  interest,  both  the  single  layer  modds 
and  the  multilayer  model  are  used.  For  the  two>layer  multilayer  model  shown  in 
Table  2,  a  third  layer  is  also  added  to  dmulate  the  actual  geometry  of  a  tree  more 
accurately  as  discussed  in  Chapter  3.  The  fidds  near  a  didectric  cylinder  using 
single  layer,  two-layer,  and  three-layer  modds  are  all  analj^ed. 

3.1  Results  for  Cylinders  of  Infinite  Length 

Most  of  the  analysis  consists  of  didectric  flinders  that  are  0.15  meters  in  radius  and 
at  a  fieequency  of  O.TSGHs.  The  didectric  constants  and  thicknesses  are  varied  for 
each  of  the  two  or  three  layers  of  the  cylinder.  The  main  trends  in  the  analysis  are 
shown  in  the  next  few  figures.  The  results  for  a  typical  three-layer  qrlinder  modd 
for  both  TM  and  TE  pohnsations  at  radial  distances  of  4,  10,  and  30  meters  are 
shown  in  Figures  16,  17,  and  18,  respectivdy.  This  three-layer  modd  has  an  inner 
layer  of  thickness  0.07m  and  ^sl5-j7,  a  second  layer  of  thickness  0.07m  and  ^=:7-j3, 
and  an  outer  layer  of  thickness  0.01m  and  er=4-jl.  These  figures  show  the  trends 
in  the  total  fidd  produced  by  a  tyjncal  tree  at  various  radial  distances.  It  is  worth 
noting  that  at  even  30  meters  (75  wavdengths)  the  total  fidd  still  shows  a  variation 
directly  behind  the.  qrlinder  about  ±1  dB  from  the  firee  space  vdue.  This  result 
is  tyincd  for  the  various  combinations  of  didectric  constants  and  thicknesses  used 
in  the  blinder  modds. 

For  comparison,  results  for  a  two-layer  cylinder  modd  at  Rs4  meters  are  shown 
in  Figure  19.  This  two-layer  modd  has  an  ixmer  layer  of  thickness  0.13m  and  6r=8-jl, 
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utd  an  outer  layer  of  thickneis  0.02m  and  er=4-j0.  The  field  plot  is  very  to 

the  previous  cate  for  the  three-layer  model  hut  differs  slightly  in  the  forward  scatter 
r^on.  This  two-layer  geometry  alto  shows  a  similar  fluctuation  in  the  field  at  30 
meters  at  the  previous  three-layer  example.  For  further  comparison,  the  dielectric 
constants  that  are  used  to  model  the  RCS  of  the  logs  in  Chapter  4  are  also  used  in 
the  near  field  analysis.  Two  tin^  layer  cylinder  geometries  are  analyised  each  with  a 
different  dielectric  constant.  A  field  plot  of  the  first  tingle  layer  model  at  R=4  meters 
is  shown  in  Figure  20  for  a  dielectric  constant  of  £r =12.4-34.9.  The  field  pattern  at 
R=4  meters  for  another  tin^e  layer  model  with  didectric  constant  £r=24-j8  is  shown 
in  Figure  21.  The  similarities  between  the  angle,  two,  and  three  layer  cylinder  models 
illustrate  the  trends  observed  in  the  fidd  patterns  of  the  various  cylinder  geometries 
analysed.  AH  the  geometries  examined  showed  approximatdy  the  same  total  fidd 
patterns  as  these  two  cases.  This  invariance  in  the  fidd  patterns  to  the  constituent 
layers  for  the  various  cases  shows  that  the  fidd  is  rdativdy  insensitive  to  the  cylinder 
layer  thickness  and  dielectric  constants  used  in  modeling  the  tree  trunk. 

To  further  illustrate  the  behaviour  of  the  fidds  with  distance,  a  rectangular  grid 
of  fidd  points  in  the  forward  scatter  rejpon  is  calculated.  Figures  22  and  23  show 
the  fidd  for  a  grid  of  points  behind  the  cylinder  for  TM  and  TE  polarisations, 
respectivdy.  The  cylinder  is  comprised  of  three-layers  with  an  inner  thiclmeaa  of 
0.08  meters  and  6r=15-j7,  a  second  layer  with  a  thickness  of  0.05  meters  and  er=6- 
j3,  and  an  outer  layer  with  a  thickness  of  0.02  meters  and  6r=4-jl.  For  comparison, 
the  fields  for  another  geometry  are  plotted  in  Figures  24  and  25.  This  cylinder  has 
an  inner  layer  thickness  of  0.08  meters  and  6r=6-j3,  a  second  layer  with  a  thirlrn^tn 
of  0.05  meters  and  6r=15-j7,  and  an  outer  la3rer  of  thickness  0.02  meters  and  6r=4- 
jl.  The  fidd  for  a  single  layer  cylinder  with  a  dielectric  constant  of  er=12.4-j4.9  is 
shown  in  Figures  26  and  27  for  TM  and  TE  polarizations,  respectivdy.  Another 
sinf^  layer  cylinder  with  a  didectric  constant  of  er=24-j8  is  plotted  in  Figures  28 
and  29,  respectivdy.  The  main  difference  in  the  plots  is  mainly  in  the  depth  of 
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Flgore  16:  Total  E-fieid  magnitude  in  dB  relative  to  the  incident  field  for  a  three- 
layer  qrliada  at  a  distance  of  4  meters  at  0.75GHs. 
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Figure  17:  Total  E-field  magnitude  in  dB  zeUtive  to  the  incident  field  for  a  three- 
layer  cylinder  at  a  distance  of  10  meters  at  0.75GHz. 
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Fignxe  18:  Total  E-field  magnitude  in  dB  relative  to  tlie  incident  field  for  a  three- 
layer  cyliiider  at  a  distance  of  30  meters  at  0.75GHz. 
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Figure  19:  Total  E-field  magnitude  in  dB  relatiye  to  the  incident  field  for  a  two-layer 
cylinder  at  a  distance  of  4  meters  at  0.75GHz. 
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Figure  20:  Total  £-fldd  magnitude  in  dB  relative  to  the  incident  field  for  a  single 
layer  cylinder  with  er=:12.4-j4.9  at  a  distance  of  4  meters  at  0.75GHz. 
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TM  polarization 


f  (degrees) 

Figure  21:  Total  E-field  magnitude  in  dB  relative  to  the  incident  field  for  a  single 
layer  cylinder  with  6r=24-j8  at  a  distance  of  4  meters  at  0.75GHz. 
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the  scattered  field  near  the  cylinder.  The  similarity  of  the  plots  for  their  respective 
polarications  again  shows  the  same  relative  insensitivity  in  near  field  patterns  to  the 
dielectric  constant  and  the  layer  thickness  as  was  observed  earlier.  And  since  the 
fields  are  relatively  insensitive  to  the  number  of  layers  used,  modeling  the  trunk  with 
a  single  dielectric  layer  is  also  valid  for  the  near-zone  scattered  field  as  is  found  for 
the  far-zone  RCS. 

4  Results  For  Cylinders  Of  Finite  Length 

The  previous  section  describes  the  results  for  infinite  length  dielectric  cylinders  un¬ 
der  normal  incidence.  The  infinite  length  of  the  cylinder  used  to  model  the  tree 
trunk  isn’t  very  faithful  to  the  true  finite  length  of  a  tree  trunk.  For  this  reason, 
the  scattered  field  of  finite  length  dielectric  cylinder  at  oblique  incidence  is  also  of 
interest.  However,  no  solution  is  readily  available  for  the  near  field  of  a  finite  length 
multilayered  dielectric  cylinders  at  oblique  incidence.  Therefore,  it  is  desirable  to 
find  a  range  over  which  the  exact  solution  for  normal  incidence  of  infinite  length 
cylinders  can  be  applied  to  the  case  of  finite  length  cylinders  with  oblique  incidence. 
To  accomplish  this,  a  UTD  solution  for  oblique  incidence  is  used  for  comparison 
at  various  inclination  angles,  d  (see  Figure  14),  with  the  exact,  normal  incidence 
solution.  However,  the  only  accurate  UTD  solution  for  the  deep  shadow  repon 
presently  available  is  for  PEC  cylinden.  The  exact  solution  for  infinite  length  cylin¬ 
ders  at  normal  incidence  is  used  to  calculate  the  field  at  various  distances  R  from 
a  PEC  cylinder.  The  UTD  solution  is  then  used  to  calculate  the  field  for  various 
incidence  angles  0  at  the  same  distances. 

For  the  geometries  analyzed,  a  cylinder  length  of  L=10  meters  and  a  radius  of 
r=:0.15  meters  as  shown  in  Figure  14  are  used.  The  two  solutions  match  very  closely 
for  the  case  of  normal  incidence  (d=0  degrees)  for  the  finite  length  cylinder.  This 
shows  that  the  dif&action  from  the  cylinder  ends  is  relatively  insignificant  for  the 
hei^t  and  radius  analyzed.  Comparison  plots  between  the  exact  solution  and  the 
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Fignxe  23:  Total  E-fidd  magnitude  in  dB  zeUtive  to  the  incident  fidd  for  a  three- 
layer  cylinder  at  O.TSGHa  for  TE  polarisation. 
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Figure  24:  Total  £-fidd  magnitude  in  dB  xdative  to  the  inddent  field  for  a  three* 
layer  cylinder  at  0.75GHz  for  TM  pdaxizadon. 
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UTD  scdution  aie  overlayed  in  Figures  30  and  31  at  a  distance  of  5  meters  for  15 
and  30  degree  indinations,  respectively.  TKe  two  solutions  are  also  overlayed  for 
a  distance  of  10  meters  and  an  inddent  angle  of  15  degrees  in  Figure  32.  The 
overlayed  plots  show  a  very  dose  similarity  in  their  near  fidd  patterns  except  for  a 
slight  variance  in  the  siddobe  spacing  and  a  difference  in  the  depth  of  the  main  null 
in  the  forward  scatter  region.  These  differences  are  relativdy  small  and  show  that  the 
normal  solution  can  be  used  for  oblique  inddence  up  to  30  degrees  depending  on  the 
accuracy  required.  However,  the  solution  for  normal  inddence  is  no  longer  applicable 
to  the  oblique  inddence  geometry  as  the  inddence  angle  B  increases  and/or  the  radial 
distance  R  increases  since  then  the  direct  ray  dominates  the  near  field  pattern. 

The  case  for  the  cylinder  above  an  infinite  ground  plane  is  analyzed  next  using 
the  same  NEC*BSG  code  as  the  above  UTD  solution.  The  case  for  a  cylinder  radius 
of  0.15m,  length  of  L=10  meters,  and  height  of  H=:2  meters  for  various  inclination 
angles  as  shown  in  Figure  15  are  analyzed.  The  near  field  patterns  for  TM  and  TE 
p<daiization8  at  a  distance  of  5  meters  for  inddence  angles  of  15  and  30  degrees 
are  overlayed  with  the  exact  solution  for  normal  inddence  in  Figures  33  and  34, 
respectively.  The  imaged  UTD  oblique  inddence  solution  and  the  normal  inddence 
exact  solution  agree  very  dosely  in  the  same  regions  and  disagree  in  the  same  re^ons 
as  the  above  exact  and  UTD  solution  without  the  ground  plane. 

5  Summary 

The  scattered  fields  in  the  vicinity  of  a  cylinder  used  to  model  a  tree  trunk  are  una¬ 
ided.  It  was  shown  that  the  scattered  fidds  for  a  infinite  length  dielectric  cylinder 
are  relatively  insensitive  to  the  dielectric  constants  and  layering  of  the  qrlinder  and 
that  a  single  layer  model  is  also  valid  for  the  near-zone  scattering.  Also,  the  blockage 
in  the  forward  scatter  caused  by  the  tree  trunk  is  shown  to  be  large.  The  blockage 
is  apprmdmately  5dB  directly  behind  the  trunk  and  decreases  to  approximatdy  IdB 
at  30  meters  from  the  trunk.  The  solution  of  the  infinite  length  didectric  cylinder 
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Figure  30:  Comparison  of  total  £-fidd  magzutude  in  dB  reUtive  to  the  incident  field 
fm  a  P.E.C.  cylinder  for  an  exact  solution  at  normal  incidence  and  a  UTD  solution 
at  15  degrees  incidence  at  a  distance  of  5  meters  at  0.75GHz 
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Figure  31:  Comparison  of  total  E-field  m^nitude  in  dB  relative  to  the  incident  field 
for  a  P.E.C.  cylinder  for  an  exact  solution  at  normal  incidence  and  a  UTD  solution 
at  30  degrees  incidence  at  a  distance  of  5  meters  at  0.75GHz. 
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Figure  32:  Comparison  of  total  £-fidd  magnitude  in  dB  relative  to  the  incident  field 
for  a  P.E.C.  cylinder  for  an  exact  solution  at  normal  inddence  and  a  UTD  solution 
at  15  degrees  inddence  at  a  distance  of  10  meters  at  0.75GHz. 
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Figure  33:  Compaiison  of  total  £-iield  magnitude  in  dB  relative  to  the  incident  field 
for  a  P.E.C.  cylinder  for  an  exact  solution  at  normal  incidence  and  an  im^ed  UTD 
solution  at  15  degrees  incidence  at  a  distance  of  5  meters  at  0.75GHz. 
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Figuie  34:  Comparison  of  total  £-firid  magnitude  in  dB  relative  to  the  incident  field 
for  a  P.E.C.  cylinder  for  an  exact  solution  at  normal  incidence  and  an  imaged  UTD 
solution  at  30  degrees  incidence  at  a  distance  of  5  meters  at  0.75GHz. 
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is  then  shown  to  be  applicable  to  a  finite  length  dielectric  cylinder  as  long  as  the 
incidence  angle  is  less  than  30  degrees  and  no  direct  ray  between  the  source  and  re¬ 
ceiver  exists.  Similar  results  are  found  for  the  finite  length  cylinder  above  a  ground 
plane  as  for  the  case  without  the  ground  plane. 
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SECTION  6 


Trunk-Ground  Interactions 

1  Introduction 

The  effect  of  the  ground  on  the  scattering  by  a  tree  trunk  is  examined.  This  builds 
on  the  results  of  Chapter  5  which  examines  the  fields  in  the  vicinity  of  infinite  length 
cylinders  and  for  finite  length  cylinders  above  a  perfect  ground  plane.  The  analy¬ 
sis  in  Chapter  5  uses  an  exact  eigenfunction  solution  and  a  UTD  analysis  with  the 
NEC-BSC  code.  In  this  chapter,  the  scattering  from  only  the  interaction  between 
the  tree  trunk  and  the  ground  is  examined.  The  scattering  from  branches  is  dis¬ 
cussed  separately  in  Chapter  7  to  simplify  the  analysis.  The  Geometric  Theory  of 
Diffraction  (GTD)  is  used  to  find  the  scattering  of  a  finite  length  dielectric  cylinder 
representing  the  trunk  protruding  from  an  infinite  dielectric  half-space  representing 
the  ground.  The  solution  is  then  verified  against  the  PEC  top-hat  solution  formu¬ 
lated  by  N.  Akhter  [25]  and  is  also  verified  using  image  theory  for  a  dielectric  cylinder 
above  a  PEC  ground  plane.  The  GTD  solution  is  then  used  to  analyze  the  scattering 
from  the  trunk-ground  interactions  for  various  trunk  and  ground  parameters.  In  the 
analysb,  the  ground  is  modeled  as  described  in  Section  2  and  the  wood  u  modeled 
based  on  the  results  of  Chapters  4  and  5. 
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Figure  35:  Geometry  used  to  model  tnmk'ground  interactions. 

2  Modeling  of  Geometry 

In  Section  2,  only  the  large  ^'anches  and  the  tnink  need  to  be  accounted  for  in 
modeling  the  tree  for  the  frequendes  between  0.3-1.3GHa.  In  Chapters  4,  5,  and  7 
it  is  valid  to  model  the  scattering  from  a  log  (wood)  as  a  dielectric  cylinder.  The 
same  method  is  used  in  this  analysis.  The  ground  is  now  modeled  as  a  dielectric  half¬ 
space  as  mentioned  in  Section  4  which  is  more  realistic  than  modeling  it  as  a  PEC 
as  is  done  in  Chapter  5.  This  modeling  of  the  geometry  will  enable  a  more  accurate 
anal3rsis  of  the  RCS  of  a  tree  by  induding  trunk-ground  interactions  (TGI).  The 
resulting  geometry  is  shown  in  Figure  35.  The  geometry  consists  of  a  finite  length 
dielectric  cylinder  representing  a  tree  trunk  protruding  from  a  didectric  half-space 
whidi  simulates  the  ground.  The  cylinder  can  consist  of  dther  multiple  concentric 
dielectric  layers  or  as  a  single  dielectric  piece.  The  polarization  defined  as  horizontal, 
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(3)  RO 
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Figure  36:  Diffiraction  interactions  used  in  trunk-ground  interaction  modeling 

with  the  electric  field  normal  to  the  x-z  plane,  and  vertical,  with  the  electric  field 
parallel  to  the  x-z  plane. 

The  modified  Geometrical  Theory  of  Diffraction  for  dielectric  materials  is  then 
used  to  calculate  the  scattering  from  this  geometry.  Due  to  rotational  symmetry, 
the  backscattered  field  needs  only  to  be  calculated  in  the  plane  containing  the  axis 
of  the  cylinder.  The  only  rays  that  need  to  be  considered  to  find  the  backscattered 
field  are  shown  in  Figure  36. 

Diffraction  from  the  cylinder-ground  junction  is  not  considered  in  this  analyris. 
This  is  because  when  the  cylinder  and  base  are  at  right  angles  and  are  of  the  same 
material,  the  scattering  from  this  junction  will  be  zero.  When  the  cylinder  and 
base  are  of  two  different  materials,  it  is  felt  that  the  scattering  from  this  interface 
will  be  small  when  compared  to  the  other  scattering  mechanisms.  Furthermore,  an 
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appropriate  Mdution  to  the  tcattering  from  this  inter&ce  is  not  presently  available. 
Therefore  without  significant  impact  on  the  calculation  of  the  backscattered  field, 
this  contribution  can  be  ignored. 

The  remaining  terms  in  the  GTD  solution  can  readily  be  calculated  using  Equa¬ 
tions  (2)  and  (12).  Accordingly,  for  the  sin^e  dificaction  Di  terms,  the  incident  and 
reflected  angles  referred  to  the  top  of  the  cylinder  are  given  by: 

r  =  ir-« 

9  =  x-e-€.  (25) 

Where  a  small  bistatic  an^e  e  is  subtracted  from  the  dif&acted  ray  to  avoid  compu¬ 
tational  errors  when  the  difEracted  ray  lies  along  the  incident  ray  path.  The  phase 
length  of  ray  Dl  referenced  to  the  center  of  the  cylinder  base  is: 

p  =  H  ^8  +  a  cos  9.  (26) 

The  argument  of  the  Ghiang  transition  function  U  is: 


U  =  2hacos6. 


The  quantity  pven  by  Equation  (5)  is: 


^  2cos6' 


Combining  Equations  (25),  (26),  (27)  ,  and  (28)  with  Equations  (2)  and  (12)  yields 
the  diffracted  field  for  location  1  given  by: 


tin  S4«  CM  S) 


2co86  *' 


I  cos  6) 


where: 


' 

E.  Boris 

' 

H,  Vcrti 


Horizontal  polarization 
Vertical  polarization. 
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Likewise  using  the  same  procedure  the  diifiracted  field  for  location  2  is  found  to 


be: 

(^^^)]) 

g2ifc(hdns+«c«.s)^-}  /IZIl7;(_2fca  cos  $).  (31) 

V  2 COSO 

For  the  Reflected-Di&acted-Reflected  (RDR)  term,  a  giTnilar  procedure  is  used 
as  before  except  the  reflection  coefficient  from  the  ground  Rg  must  be  taken  into 
account.  The  resulting  diffracted  fidd  for  the  RDR  term  is: 

w)’ { (“•  (^)  *  -  ("^)|  * 

e^c*^*^“**y^^^Te(2fcocos^).  (32) 

The  Diffracted-Reflected  (DR)  and  Reflected-Di&acted  (RD)  terms  can  be  cal¬ 
culated  as  follows.  It  has  been  shown  that  the  DR  and  RD  terms  give  the  complete 
solution  in  the  dihedral  region  if  a  small  bistatic  angle  is  taken  between  the  inddent 
and  reflected  ray  paths  for  both  the  RD  and  DR  terms  [25].  The  appropriate  ray 
paths  are  shown  in  Figure  37.  When  the  RD  and  DR  terms  are  calculated  this 
way  and  combined,  the  discontinuities  in  thdr  solutions  will  cancel  and  provide  the 
correct  bounded  result.  This  combination  is  done  numerically  by  the  computer  and 
each  dii&action  term  is  presented  independently.  Using  the  geometry  of  Figure  37, 
the  DR  field  is: 


{ [“•  (^) + (^)]  + 

gi*(a co<S+«)-A.m(fl+e))g^^+»*  eo.(«-tan  ‘  /_£_T(2fca  COS  ^) 

V  2cose  '  ' 
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and  the  RD  field  is 


re(2Jfeacos«).  (34) 

V  2  cos  9 

It  should  be  noted  that  the  method  to  calculate  the  diffracted  fields  from  the  TGI 
geometry  doesn’t  take  the  creeping  wave  component  of  the  cylinder  into  account. 
The  above  di&acted  fields  are  combined  to  calculate  the  scattered  field  for  the 
trunk-ground  interactions. 

3  Verification  of  TGI  Solution 

The  GTD  solution  for  the  trunk-ground  interaction  is  verified  in  two  ways.  The  RCS 
of  a  finite  length  dielectric  cylinder  is  calculated  using  a  modified  exact  eigenfunction 
expansion.  Then  the  RCS  of  a  dielectric  cylinder  of  half  the  length  as  used  in  the 
exact  solution  above  an  infinite  ground  plane  is  analyzed  using  the  TGI  solution  for 
9  approximately  0  degrees.  The  two  values  agree  quite  well  for  the  various  cases 
tried. 

Unfortunately,  no  other  solution  for  the  dielectric  trunk-ground  interaction  is 
available,  but  a  solution  is  for  a  PEC  top-hat.  The  top-hat  has  a  finite  radius 
PEC  disk  for  a  base  instead  of  an  infinite  ground  plane  as  in  the  TGI  geometry. 
Therefore,  to  compare  the  two  solutions,  the  equivalent  top-hat  geometry  must  be 
chosen  carefully.  First,  the  diffraction  from  the  edge  of  the  base  of  the  top-hat 
must  be  eliminated  since  this  edge  doesn’t  exist  for  the  TGI  geometry.  Second, 
the  dimensions  of  the  top-hat  must  be  chosen  such  that  only  RD,  DR,  and  RDR 
diffractions  can  occur  from  the  top  of  the  top-hat  cylinder  and  not  from  the  base 
edge.  Once  these  modifications  are  done,  the  solutions  can  be  verified  for  the  PEC 
case.  This  is  shown  in  Figure  38  for  a  top-hat  of  radius  a=0.2m  and  height  H=15m. 
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Tkble  4:  Diameters  of  trees  {DBH)  and  the  corresponding  tree  height  h  and  effective 
tree  height  ht  of  the  cylinder  used  to  model  the  tree. 


DBH 

tree  height  (h) 

effective  tree  height  {ht) 

12”  (0.3048m) 

17.4m 

14.8m 

14”  (0.3556m) 

20.4m 

17m 

16”  (0.4064m) 

23.3m 

20m 

The  base  radius  of  the  top-hat  is  taken  as  430  meters  to  eliminate  the  RD,DR,  and 
RDR  diffraction  terms  from  occurring  from  the  base  edge.  Since  a  solution  for  a 
dielectric  top-hat  is  not  readily  available,  the  verification  of  the  solution  is  complete 
and  the  solution  of  the  TGI  geometry  is  assumed  correct  for  the  dielectric  case  based 
on  the  above  two  cases. 

4  Analysis  of  Trunk- Ground  Interactions 

The  trunk-ground  interactions  are  analyzed  unng  the  above  GTD  TGI  solution  for 
the  scattered  fields.  For  the  tree  geometries,  the  diameters  of  the  trees  are  taken 
from  measurements  at  a  forested  area  near  OSU’s  Big  Ear  telescope.  The  trees  in 
the  forest  range  in  diameters  from  small  saplings  up  to  trees  16"  in  diameter.  The 
analysis  conasts  mostly  of  trees  with  DBHh  between  These  dimensions 

are  large  enough  for  a  GTD  analysis  and  comprise  a  fairly  large  number  of  the  trees 
in  the  forest.  Using  these  diameters  and  Equations  (14)  and  (17)  appropriate  heights 
of  the  trees  can  be  found.  The  dielectric  constants  are  chosen  based  on  the  results  of 
Chapter  4  for  the  tree  and  on  Section  2  for  the  ground.  The  results  of  the  analysis 
axe  described  below. 

The  RCS  for  a  tree  with  DBH=li”  and  a  corresponding  hdght  from  Table  4 
of  17  meters  is  shown  in  Figure  39  at  a  frequency  of  0.75GHz.  In  the  Figure,  the 
top  plot  is  for  a  cylinder  with  a  dielectric  constant  of  6(=12.4-j4.9  while  the  bottom 
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Fignxe  38:  Comparison  of  top-hat  and  TGI  solutions  foi  the  PEC  case.  The  cylin¬ 
der  radins  a=0.2m  and  height  H=15m.  The  top-hat  base  is  430  meters  for  this 
comparison  and  the  frequency  is  0.75GHz. 


plot  is  for  a  cylmdei  with  a  dielectric  constant  of  et=24.8-j8.  The  ground  for  both 
plots  has  a  dielectric  constant  of  ey=20-j0.15  which  from  Table  3  corresponds  to 
the  measured  soil  at  Ohio  State  University’s  Big  Ear  telescope.  The  same  size 
tree  is  shown  in  Figure  40  and  Figure  41  for  soil  dielectric  constants  of  ep=7-jl.5 
and  6^=18*35,  respectively.  The  effect  of  the  soil  and  trunk  dielectric  constants  sire 
apparent  in  the  PEC  case  of  Figure  38  where  a  relatively  flat  RCS  with  incidence 
angle  occurs.  The  dielectric  cases  show  a  large  null  in  the  pattern  around  0=15 
degrees  and  another  at  around  0=75  degrees  for  only  the  vertical  polarization  but 
no  nulls  for  the  horizontal  polsirization.  These  angles  correspond  to  the  Brewster 
angle  for  the  verticsd  polarization  for  the  dielectric  surfaces  of  the  trunk  smd  ground. 
No  nulls  are  present  for  the  horizontal  polarization  since  the  Brewster  angle  for 
the  horizontal  polarization  is  dependent  on  the  magnetic  properties  of  the  material. 
Since  the  trunk  and  ground  are  assumed  to  have  the  permeability  of  free  space  (/io) 
the  Brewster  angle  doesn’t  exist  for  the  horizontal  polarization. 

The  difference  in  the  polarizations  is  worth  further  examination.  The  reflection 
coeffidents  from  a  dielectric  half-space  are  shown  in  Figure  42  for  the  dielectric 
constants  used  to  model  the  ground  and  wood.  The  reflection  coeflidents  for  the 
ground  (F^)  show  a  deep  null  around  15  degrees.  The  location  and  depth  of  the  null 
depend  on  the  exact  dielectric  material  of  interest.  A  mmilar  null  is  observed  in  the 
reflection  coeffidents  for  the  tree  dielecti  materials  (Ft).  The  reflection  coefiident 
for  the  tree  dielectric  constants  is  shown  for  a  dielectric  half-space  instead  of  for 
a  cylinder.  This  is  for  ease  of  computation  and  to  show  the  inherent  behaviour  of 
the  dielectric  material  so  the  effect  of  the  cylindrical  geometry  is  removed.  The 
inaccuracy  caused  by  this  is  small  since  the  material  is  lossy  so  any  contribution  due 
to  reflections  from  the  back  face  of  the  cylinder  would  be  attenuated  and  the  tnain 
contribution  will  come  from  the  front  cylinder  face  which  is  BitnilaT  to  the  case  for 
a  half-space.  The  reflection  coeffident  for  the  vertical  polarization  shows  a  BimilaT 
null  for  the  tree  dielectric  materials  compared  to  the  soil  dielectric  materials  but 
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Figure  39:  Radar  cross  section  for  a  tree  with  DBH=14”  and  height  ht=l7m.  at  a 
frequency  of  0.75GHz  for  both  horizontal  and  vertical  polarizations.  The  relative 
dielectric  constant  of  the  ground  is  ^=20-j0.15.  The  rdative  dielectric  constant  of 
the  tree  is  e(=12.4-j4.9  for  the  top  plot  and  ct=24-j8  for  the  bottom  plot. 
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Fignxe  40:  Radar  cross  section  for  a  tree  with  DBH=li”  and  height  ht=17m  at  a 
frequency  of  0.75GHz  for  both  horizontal  and  vertical  polarizations.  The  relative 
dielectric  constant  of  the  ground  is  ^=7-jl.5.  The  relative  dielectric  constant  of  the 
tree  is  e(=12.4-j4.9  for  the  top  plot  and  e(=24>j8  for  the  bottom  plot. 
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Figure  41:  Radar  cross  section  for  a  tree  with  DBH=14”  and  height  ht=17m  at  a 
frequency  of  0.75GHz  for  both  horizontal  and  vertical  polarizations.  The  relative 
dielectric  constant  of  the  ground  is  ^=18*j5.  The  relative  dielectric  constant  of  the 
tree  is  e(=12.4-j4.9  for  the  top  plot  and  e(=24-j8  for  the  bottom  plot. 
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at  higher  incident  angles  from  horiaontal  since  the  dielectric  half-space  for  the  tree 
dielectric  constants  is  oriented  vertically. 

These  nulls  in  the  reflection  coefficient  significantly  affect  the  RCS  of  the  tree 
geometry  since  the  dominant  contribution  to  the  RCS  of  the  TGI  geometry  is  from 
the  RD  and  DR  terms  which  represent  the  dihedral  effect  between  the  trunk  and 
the  ground.  Since  the  RD  and  DR  terms  both  are  multiplied  by  the  reflection  from 
the  ground,  a  decrease  in  the  ground  reflection  coefficient  will  significantly  decrease 
the  RD  and  DR  term  magnitudes  near  the  ground  Brewster  angle.  This  decrease  in 
the  reflection  coefficient  is  only  for  the  vertical  polarization  and  not  the  horizontal 
polarization.  The  same  effect  occurs  for  the  RDR  term  but  the  effect  of  the  RDR 
term  is  smaller  than  the  effect  of  the  RD  and  DR  terms.  As  the  incidence  angle 
increases,  the  Brewster  angle  from  the  tree  surface  has  a  significant  effect  on  the 
dif&acted  fields.  This  can  be  seen  by  examining  Equation  (12)  and  noting  that  the 
reflection  coeffident  from  the  dielectric  surface  (the  trunk)  modifies  the  diffraction 
coeffidents.  Therefore,  as  the  tree  reflection  coeffident  decreases,  the  diffracted  field 
will  decrease  which  will  affect  all  the  di&action  terms.  Therefore,  we  can  condude 
that  the  effect  of  the  reflection  coeffidents  near  the  Brewster  angles  on  the  reflected 
fidd  from  the  didectric  surfaces  causes  a  significant  difference  between  the  RCS  of 
two  polazizations  from  the  tree. 

To  further  illustrate  the  effect  of  the  didectric  materials  on  the  RCS  of  the 
tree,  other  cases  will  be  examined.  In  Figure  43,  the  RCS  at  0.75GHz  for  a  tree 
of  DBH=12"  with  a  corresponding  effective  height  ht  from  Equation  (14)  of  11.6 
meters.  The  didectric  constants  of  the  tree  are  et=12.4-j4.9  for  one  graph  and 
e(=24.8-j8  for  the  other  graph.  The  ground  didectric  constant  is  ^=18-j5  for  both 
plots.  The  RCS  shows  the  same  pattern  as  shown  before  for  trees  with  diameters 
of  16".  This  shows  that  the  diameter  has  Httle  effect  on  the  difference  between 
polazizations.  In  Figure  44  the  same  size  tree  and  same  dielectric  constants  as 
in  Figure  43  aie  used  but  this  time  at  a  frequency  of  1.3GHz  which  is  the  highest 
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Figue  42:  Reflection  coeflident  from  didectzic  half-spaces  to  show  the  difierences 
in  the  reflection  coefiBdent  (F)  for  the  horisontal  and  vertical  polarisations  of  the 
dielectric  materials  need  in  modehng  the  tmnk-gronnd  interaction.  Top  plot  is  of 
the  reflection  coefficients  for  sdl  didectric  constants  (F^)  and  the  bottom  plot  is  of 
the  reflection  coeffidents  for  tree  didectric  constants  (Ft). 


frequency  in  the  band  of  interest.  Figure  45  shows  the  RCS  for  a  tree  with  DBH=16” 
and  a  height  of  20in  at  a  frequency  of  1.3GHz.  The  dielectric  constant  of  the  ground 
is  60=7*jl.5.  The  same  difference  between  the  horusontal  and  vertical  polarizations 
is  seen  for  these  tree  geometries  as  in  the  previous  cases. 

5  Summary 

The  analysis  of  the  scattering  from  a  tree  trunk  including  the  effect  of  the  ground 
is  examined.  The  trunk  is  modeled  as  a  finite  length  dielectric  cylinder  and  the 
ground  is  modeled  as  an  infinite  dielectric  half-space.  The  effect  of  the  dielectric 
constants  for  the  ground  and  tree  is  examined.  It  is  found  that  for  all  the  dielectric 
constants  chosen,  the  RCS  of  the  tree  for  horizontal  polarization  is  5*7dB  higher 
than  for  the  vertical  polarization.  This  is  explained  by  the  Brewster  angle  of  the 
dielectric  materials  (both  the  tree  and  ground).  At  the  Brewster  angle,  the  reflection 
from  the  dielectric  surfrces  is  much  higher  for  the  horizontally  polarized  field  than 
the  vertically  field.  The  reflection  coefSdent  for  the  vertical  polarization  is  close  to 
zero  at  the  Brewster  angle.  This  difference  in  reflection  coefiicients  for  the  polariza¬ 
tions  causes  the  large  difference  in  the  RCS  of  a  tree  for  the  horizontal  and  vertical 
polarizatioru. 
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9  (decrees) 

Fignze  43:  Radar  ciom  lection  for  a  tree  with  DBH  =  I2f'  and  height  h«=14.8m  at 
a  frequency  of  0.75GHa  for  both  honzontal  and  vertical  polariaatiouB.  The  rdative 
dielectric  conitant  of  the  ground  is  18-j5.  The  rdative  didectric  constant  of  the 
tree  is  Ct=12.4-j4.9  for  the  top  plot  and  6t=:24-j8  for  the  bottom  plot. 
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9  (defrees) 

Fignze  44:  Radar  croii  section  for  a  tree  with  DBH=12”  and  height  ht=14.8m  at 
a  frequency  of  1.3GHs  ibi  both  hozisontal  and  Tertical  polarisations.  The  relative 
dklectiic  OMutant  of  the  ground  is  ^=18-j5.  The  relative  dielectric  constant  of  the 
tree  is  e(=:12.4-j4.9  for  the  top  plot  and  6t=24-j8  for  the  bottom  plot. 
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Figure  45:  Radar  cross  section  for  a  tree  with  DBH=1&*  and  height  ht=20in  at 
a  frequency  of  OGHs  for  both  horisontal  and  vertical  polarisations.  The  relative 
dielectric  constant  of  the  ground  is  e,— 7-jl.5.  The  relative  dielectric  constant  of  the 
tree  is  e(=12.4-j4.9  for  the  top  plot  and  6t=24-j8  for  the  bottom  plot. 


SECTION  7 
Branch  Scattering 

1  Introduction 

This  chapter  discusses  the  scattering  from  tree  branches.  In  Chapter  6  the  scatter¬ 
ing  due  to  the  dihedral  between  the  trunk  and  ground  is  studied.  The  polarization 
difference  in  the  ROS  of  the  trunk-ground  is  also  shown  in  Chapter  6.  The  question 
now  is  where  does  the  level  of  scattering  &om  branches  compare  to  the  level  of  scat¬ 
tering  from  the  trunk-ground  interaction.  Due  to  the  complexity  of  tree  branches,  a 
solution  involving  all  the  scattering  interactions  would  be  very  complex.  Therefore, 
a  simplified  branch  geometry  that  fits  the  needs  of  the  analysis  must  be  decided 
upon.  The  geometry  coxisists  only  of  the  interactions  between  a  single  branch  and 
the  trunk.  This  involves  analyzing  the  scattering  from  the  dihedral  between  a  tree 
branch  and  trunk,  and  the  specular  scattering  from  the  side  of  a  tree  branch.  To  find 
the  level  of  scattering  from  branches,  it  is  desired  to  find  the  wiM-TimnTn  scattering 
from  a  single  branch  that  could  be  expected.  To  accomplish  this,  a  series  of  worst 
case  scenarios  are  analyzed  that  will  provide  the  highest  level  of  scattering  from  tree 
branches  that  would  occur.  Based  on  the  results  for  a  single  branch,  the  importance 
of  the  scattering  from  branches  can  be  deduced. 

The  analyns  consists  of  modeling  the  tree  trunk  and  branches  with  dielectric 
cylinders  as  is  done  in  Chapters  4  and  5.  The  dielectric  constants  are  chosen  based 
on  information  in  the  literature  that  is  presented  in  Table  1  and  on  the  results  of 
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Chapter  4.  The  wood  models  are  then  used  in  a  UTD  solution  for  dihedrals  of 
dielectric  cylinders.  Based  on  this  analysis,  a  level  for  the  RCS  from  a  branch  is 
found. 


2  Modeling  of  Branch  Scattering 

The  actual  geometry  of  tree  branches  is  very  complex  and  is  beyond  the  level  of 
analysis  desired.  The  purpose  of  the  analysis  is  to  find  a  level  for  the  scattering 
from  a  single  branch  and,  therefore,  dedde  on  the  importance  of  branch  scattering 
compared  to  the  trunk-ground  interactions.  Since  the  exact  level  of  branch  scattering 
is  not  of  interest,  an  attempt  will  be  made  to  simplify  the  analysis.  As  already 
discussed  in  Section  2,  the  geometry  of  the  tree  can  be  simplified  for  the  frequencies 
of  interest  to  eliminate  many  of  the  smaller  tree  branches  leaving  only  the  larger 
branches  and  the  trunk.  While  a  more  complete  modeling  of  the  branches  would 
account  for  the  various  interactions  between  the  trunk  and  the  numerous  branches, 
this  analysis  is  currently  too  complex  for  the  results  desired. 

Since  the  exact  level  of  branch  scattering  is  not  desired  and  only  the  maximum 
level  of  branch  scattering,  a  worst  case  scattering  analysis  isr  done.  To  do  this,  it  is 
observed  that  for  a  field  incident  on  the  tree,  many  of  the  larger  tree  branches  will  be 
oriented  at  obHqae  an^es  to  the  incident  field.  The  backscattered  field  for  oblique 
incidence  will  be  small  and,  therefore,  this  case  will  be  ignored.  The  dominant 
effect  will  be  when  the  incident  field  is  near  normal  incidence  (specular)  on  the 
nde  of  a  branch.  This  will  have  a  much  larger  backscattered  field  than  the  oblique 
incidence  case.  It  is  also  observed  that  the  larger  branches  of  a  tree  are  oriented  at 
many  different  angles  with  respect  to  the  vertical  axis  of  the  tree  trunk.  The  angle 
between  the  trunk  and  the  branch  form  a  dihedral  which  can  have  a  significant 
backscattered  field  if  the  incident  field  is  in  the  plane  of  the  dihedral.  The  dihedral 
effect  will  be  largest  for  backscatter  if  the  angle  between  the  branch  and  trunk  is  90 
degrees.  Using  this  reasoning,  the  only  branch  scattering  effects  of  interest  are  the 
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Figoie  46:  Geometry  used  to  model  branch  scattering. 


speculai  scattering  from  the  side  of  a  tree  branch  and  the  dihedral  effect  between  the 
trank  and  branch.  AD  other  scattering  effects  will  be  ignored  inclnding  higher  order 
interactions  between  the  branches  because  these  effects  will  be  small  compared  the 
two  main  effects  mentioned  above. 

The  geometry  that  will  be  analyzed  is  shown  in  Figure  46.  It  counts  of  two 
finite  length  dielectric  cylinders  at  a  90  degree  dihedral  angle.  For  this  geometry, 
a  horizontally  polarized  wave  is  oriented  normal  to  the  x-z  plane  while  a  vertically 
polarized  wave  is  parallel  to  the  x-z  plane.  The  backscattered  field  is  calculated 
using  the  UTD  solution  formulated  in  [25]  and  briefly  summarized  in  Appendix  C. 
Care  must  be  taken  when  uring  this  solution  to  guarantee  that  the  diameters  of  the 
dielectric  cylinders  are  greater  than  approximately  0.7  wavelengths  otherwise  the 
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UTD  solution  will  break  down.  For  this  geometry,  the  specular  scattering  from  the 
side  of  a  branch  will  occur  for  d=90  degrees. 

To  model  the  scattering  from  a  tree  branch,  the  sizes  and  dielectric  constants 
of  the  cylinders  must  be  decided  upon.  The  dielectric  constants  of  the  wood  are 
modeled  using  the  same  values  that  are  found  to  compare  with  measurements  in 
Chapter  4.  The  lengths  of  the  dielectric  cylinders  are  taken  from  Equation  (17). 
Unfortunately,  these  equations  only  provide  an  effective  tree  height  ht  and  height 
to  maximum  crown  diameter  hu.  However,  these  equations  can  be  iised  to  derive 
approximate  values  for  the  lengths  of  the  dielectric  cylinders  used  in  the  modeling. 
And  since  only  a  worst  cs^e  analysis  is  of  interest,  the  approximations  of  the  length 
are  appropriate  if  chosen  correctly. 

Therefore,  using  Equation  (17),  if  the  diameter  of  the  tree  (DBH)  is  known  then 
the  height  to  the  widest  part  of  the  canopy  hu  and  the  effective  height  of  the  tree  hf 
can  be  found.  It  is  then  assumed  that  the  longest  branches  occur  at  the  widest  part 
of  the  canopy.  While  this  assumption  may  not  always  be  valid,  we  are  assuming 
a  90  degree  dihedral  which  forces  the  longest  branches  to  occur  at  the  widest  part 
of  the  canopy.  The  length  of  the  tree  cylinder  Lt  is  then  taken  as  the  difference 
between  the  effective  height  of  the  tree  ht  and  the  widest  part  of  the  canopy  hu. 
The  diameter  of  the  cylinder  used  to  modd  the  trunk  is  taken  as  the  DBH  of  the 
tree.  These  assumptions  will  significantly  overestimate  the  actual  dimensions  of  the 
trunk  since  in  the  modd,  the  cylinder  has  a  constant  diameter  but  in  an  actual  tree 
the  trunk  will  be  tapered. 

For  the  diameter  and  length  of  the  branch  cylinder,  other  approximations  are 
used.  The  literature  is  short  on  details  on  the  modeling  tree  branches  and  formu¬ 
las  rdating  tree  idze  to  branch  sizes.  Therefore,  since  only  a  worst  case  andysis 
is  of  interest,  an  intuitive  modd  is  used  that  will  still  provide  physically  realistic 
dimensions  but  will  overestimate  the  actual  dimennons  of  the  branches  to  provide 
a  maximum  levd  of  the  RCS  &om  a  branch.  Accordingly,  it  is  assumed  that  the 
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branch  and  tree  have  approximately  the  same  dielectric  constants  which  should  be 
a  reasonable  assumption  for  the  larger  branches  which  are  the  only  ones  of  interest 
here.  The  diameter  of  the  branch  will  be  taken  to  be  equal  to  the  diameter  of  the 
trunk.  This  is  justified  rince  it  would  be  unrealistic  for  a  branch  to  have  a  larger 
diameter  than  the  tree  and  setting  the  diameters  equal  will  provide  a  worst  case 
scenario.  Choosing  an  appropriate  length  for  the  branch  is  more  difficult.  Setting 
the  branch  length  equal  to  the  trunk  length  would  provide  a  worst  case  analysis  but 
would  be  highly  unrealistic  and  nonphysical.  Therefore,  the  length  of  the  branch 
cylinder  will  be  arbitrarily  taken  as  one-half  the  length  of  the  trunk.  This  length  is 
more  realistic  than  m^lriTig  the  lengths  equal  but  should  still  provide  an  overestima- 
ii<m  of  the  RCS  of  the  branch.  Also,  the  branch  is  modeled  as  a  constant  radius, 
straii^t  cylinder  but  since  an  actual  branch  is  tapered  and  usually  bent,  the  actual 
RCS  of  a  branch  is  overestimated.  These  dimensions  for  the  branch  overestimate 
the  actual  dimensions  and,  therefore,  the  scattering  for  a  tree  branch. 

3  Analysis  of  Branch  Scattering 

The  diameters  used  for  the  tree  trunk  are  based  on  a  survey  of  trees  in  the  woods 
surrounding  the  Ohio  State  University  Radio  Telescope  known  as  Big  Ear.  The  tree 
surrqr  showed  that  there  axe  a  large  number  of  trees  with  DBH%  between  12  and 
16  inches  (0.30  and  0.41  meters).  Therefore,  the  cylinders  used  to  model  the  trunks 
and  branches  will  be  between  0.3  and  0.41  meters  in  diameter.  Using  these  values 
in  Equation  (17)  the  lengths  Lt  and  Xs  can  be  found  as  discussed  previously.  These 
dimensions  axe  summaxixed  in  Table  5  for  various  DBH  values. 

The  dimensions  from  Table  5  axe  thmi  used  with  the  dielectric  constants  from 
Chapter  4  to  find  a  level  for  the  scattering  from  tree  brandies.  The  RCS  from  the 
branch  dihedral  for  a  tree  with  a  DBHs:12”  and  didectric  constant  of  =  et=12.4- 
j4.9  is  shown  in  Figure  48  for  both  horizontal  and  vertical  polarisations  at  a  frequency 
of  0.75GHz.  For  comparison,  the  RCS  for  the  same  nze  tree  but  with  a  didectric 
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IkUe  5:  Diameters  of  trees  (DBH)  and  the  corresponding  cylinder  lengths  used  in 
modeling  branch  scattering  for  worst  case  analysis. 


DBH 

tree  cylinder  length  (X() 

branch  length  (Lb) 

12”  (0.3048m) 

5.26m 

2.63m 

14”  (0.3556m) 

6.06m 

3.03m 

16”  (0.4064m) 

6.86m 

3.43m 

constant  of  cs  s=  ei=24-j8  is  shown  in  Figure  49  at  the  same  frequency.  For  both 
Figures  48  and  49,  the  RCS  for  horisontal  polarisation  is  greater  than  the  RCS 
for  vertical  polarisation.  For  compatiscm,  a  branch  dihedral  of  the  same  sise  and 
at  the  same  frequency  is  shown  in  Figure  47  but  for  a  PEC  case.  The  PEC  case 
doesn’t  show  the  same  difference  in  the  horisontal  and  vertical  polarisations  that 
the  didectric  cases  show.  This  difference  in  the  polarisations  is  due  to  the  reflection 
coeffidents  from  the  branch  surface  being  smaller  for  vertical  polarisation  than  for 
hmisontal  polarisation.  This  is  due  to  the  decrease  in  the  reflection  coeffident  for 
vertical  polarisation  near  the  Brewster  an^e  as  is  described  in  Chapter  6.  Therefore, 
the  differences  in  the  two  polarisations  are  due  to  the  same  reasons  as  is  explained 
for  the  polarisation  differences  of  the  GTI  geometry  in  Chapter  6. 

The  RCS  of  the  branch  dihedral  should  be  examined  more  dosely.  For  0  near 
0  degrees,  most  of  the  scattering  is  from  the  trunk  cylinder  and  not  the  branch. 
The  trunk  scattering  is  discussed  in  detail  in  Chapter  6  and  isn’t  important  for  the 
current  anal^nts  of  brandi  scattering.  Near  9^90  degrees  the  scattering  is  mainly 
from  the  specular  scattering  from  the  side  of  the  branch  cylinder.  The  RCS  for  the 
branch  is  less  than  7  dBsm  for  the  geometries  in  Figures  48  and  49.  For  the  dihedral 
reipon,  the  scattering  for  the  horisontal  and  vertical  polarisations  are  approximately 
5  dB  different. 
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@  (degrees) 

Figue  47:  Radar  cron  aectkm  of  90  dq^ree  PEC  qrlmdxical  dihedral  with  diameters 
of  12”  fm  cmnparison  to  the  didectxic  can  of  scattering  from  a  tree>brandi  dihednd 
at  a  beqaency  of  0.75GHa. 
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©  (degrees) 

Figue  48:  Radar  crcwa  aection  of  90  degree  cylmdxical  dihedral  to  nmalate  the 
•catteiiiig  ficom  a  tree-brandi  dihedral  at  a  frequency  of  0.75GHb.  The  tree  diameter 
{DBS)  if  12”  with  a  relative  didectric  constant  of  =:  6(=12.4-j4.9. 
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Figue  49:  Radac  cro«t  section  of  90  d^pree  cylindrical  dihedral  to  rimnlate  the 
scattering  {coin  a  ttee>braach  dihedral  at  a  frequency  of  O.TSGHa.  The  tree  diameter 
{DBS)  is  12”  and  a  relative  dielectric  constant  of  =  Cf— 24-j8. 
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To  farther  illuttrete  the  scattering  from  branches,  several  more  geometries  are  an- 
aljrsed.  The  RCS  for  a  tree  with  didectric  constants  of  =  e(=12.4-j4.9 

and  64  s  6ts24-j8  is  shown  in  Figare  50  for  both  horisontal  and  vertical  polarisa¬ 
tions  at  a  frequency  of  0.75GHs.  In  Figure  51,  the  RCS  at  0.75GH2  of  a  tree  with 
DBJfslO”  and  with  dielectric  constants  of  es  =  e(=:l2.4-j4.9  and  cs  =  et=24-j8  is 
shown.  In  Figure  52,  the  RCS  of  the  same  tree  with  DBH=16”  and  with  didectric 
constants  of  cs  =  C|=sl2.4-j4.9  and  cs  =  et=24-j8  is  shown  but  at  a  frequency  of 
1.3GHs  which  is  the  highest  freqnenqr  in  the  band  of  interest.  By  comparing  Fig¬ 
ures  48,  49,  50,  and  51  it  is  seen  that  the  sise  of  the  branch  (diameter  and  length) 
doesn*t  affect  the  shape  of  the  RCS  considerably.  The  use  mainly  causes  the  RCS 
of  the  branches  to  shift  upward  for  larger  sise  branches.  The  difference  between 
the  RCS  for  the  horisontd  and  vertical  pdarisations  is  nearly  independent  for  the 
branch  siaes  analysed  and  is  always  approadmatdy  5  dB.  The  largest  RCS  in  the 
dihedral  region  of  the  geometnes  analysed  is  for  the  geometry  of  Figure  52.  The 
RCS  is  less  than  approximatdy  5dBsm  in  the  dihedral  refpon  and  is  less  than  12 
dBam  for  the  specular  scattering  from  the  log.  The  log  for  Figure  52  is  16**  in  di¬ 
ameter  which  is  the  largest  tree  sise  in  the  woods  at  OSU’s  radio  telescope  and  the 
frequency  is  the  highest  in  the  band  of  interest.  Since  the  dimensions  used  in  the 
analyris  are  the  largest  that  are  reasonable  iot  the  worst  case  anafyris,  the  case  in 
Figure  52  is  representative  of  the  largest  RCS  for  branch  scattering  for  the  tree  sises 
and  frequendes  of  interest. 

4  Summary 

A  model  for  the  scattering  from  the  dihedral  framed  by  the  tree  trunk  and  a  branch 
is  presented.  The  anafysis  consists  of  a  worst  case  analysis  to  show  the  highest  level 
of  brandi  scattering  that  could  be  expected.  As  a  result,  the  sixes  of  the  cylinders 
used  to  represent  the  branches  and  trunk  are  {ncked  as  large  as  conndered  msonable 
based  on  fromnlas  fra  tree  heigiht  and  on  intuition  about  branch  uses.  The  resulting 
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RCS  (dBsm)  RCS  (dBsm) 


FSgue  50:  Radar  croat  section  of  90  degree  cylindrical  dihedral  to  simnlnte  the 
scaltetiBg  from  a  tree>brancli  dihedral  at  a  frequency  of  0.75GHZ.  The  tree  diameter 
(DBS)  is  14”  with  relatiTe  dielectric  constants  of  ~  ^=12.4-j4.9  and  =  6fS=24- 


0  (degrees) 


Figuse  51:  Rader  croM  lectkm  of  90  degree  cyHsdrical  dihedral  to  rimnlate  the 
scattetmg  fram  a  tree*bra]ich  dihedral  at  a  frequency  of  0.75GEs.  The  tree  diameter 
{DBS)  u  16”  with  rdatiwe  dielectric  coxutaats  of  Cb  ===  e(sl2.4-j4.9  and  =  et=24- 
j8. 
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0  (degrees) 

Hguze  52:  Radar  cross  sectkm  of  90  degree  cyliadrical  dihedral  to  simiilate  the 
scatteriag  a  tree>hraach  dihedral  at  a  frequency  of  1.3GHs.  The  tree  diameter 
(J?BE)  is  16”  with  relatiTe  dielectric  constants  of  Cb  =  et=12.4-j4.9  and  et  =  €(=24- 

j«. 


model  ov«refiimated  the  RCS  of  actual  branches  since  in  modeling,  the  branches 
and  trunk  are  modeled  as  perfectly  straight  and  of  constant  diameter  whereas  for 
actual  trees  the  branches  and  trunk  are  bent  and  tapered  so  the  actual  geometry  will 
have  a  lower  RCS  than  modeled.  Using  this  worst  case  analysis,  the  maximum  RCS 
for  a  tree  branch  is  less  than  5  dBsm  in  the  dihedral  repon  between  the  trunk  and 
branch  and  would  be  much  leas  out  of  the  dihedral  plane.  The  specular  scattering 
from  a  branch  is  less  than  12  dBsm.  While  the  specular  scattering  from  a  branch  is 
greater  than  the  dihedral  r^on,  the  chances  of  the  incident  field  impinjpng  exactly 
perpendicular  to  a  branch  is  small  and,  therefore,  the  scattering  from  a  tree  branch 
win  be  smaller  than  calculated.  Therefore,  it  is  shown  that  the  highest  level  of 
scattering  from  branches  is  approximately  20dB  lower  than  the  scattering  from  the 
trunk-ground  interaction  in  the  dihedral  r^on. 
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SECTION  8 

Summary  and  Conclusions 


The  ekctromagnetic  scattering  from  the  trunk  and  branches  of  a  tree  is  analyzed. 
To  scdve  for  the  scattered  field,  the  modified  Geometrical  Theory  of  Oif&action  is 
used.  Befne  the  scattering  from  tree  can  be  analyzed,  some  background  information 
iz  required.  A  summary  of  the  current  techniques  and  methods  of  modeling  trees  is 
presented.  The  geometry  used  for  a  low-frequeacy  model  of  the  trunk  and  branches 
is  shown.  The  corresponding  dielectric  information  for  modeling  wood  is  presented. 
But  before  the  trunks  and  branches  are  analyzed,  the  scattering  from  tree  logs  is 
examined.  These  measurements  are  used  to  find  a  set  of  values  that  are  presented  in 
the  summary  that  can  appropriately  be  used  to  model  the  types  of  wood  of  interest. 
A  set  of  didectric  constants  for  wood  are  found  and  used  in  modeling  logs  to  dosdy 
match  the  measured  RCS.  It  is  also  found  that  the  sur&ce  roughness  can  be  ignored 
for  firequendes  bdow  apprommatdy  6GHz  for  the  trees  of  interest.  A  method  is 
presented  that  could  be  used  to  match  the  measured  RCS  of  a  log  to  theory  using 
a  roughness  modd  for  the  wood  surface.  But  due  to  sensitivity  in  the  modd  to  the 
sur&ce  roughness  and  the  variation  of  the  RCS  of  the  same  section  of  a  trees  trunk, 
it  is  shown  that  this  method  isn’t  recommended  in  modding  the  RCS  from  a  tree 
trunk  above  6GHz. 

The  results  for  the  modeling  of  a  tree  are  combined  with  the  information  pre¬ 
sented  to  calcdate  the  scattered  fidds  in  the  ddnity  of  a  tree  trunk.  It  is  shown 
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that  the  scattered  field  near  a  tree  trunk  is  relatively  insensitive  to  the  constituent 
parameters  oi  the  tree.  Also,  the  blockage  in  the  forward  scatter  re|pon  caused  by 
the  trunk  is  shown  to  be  large  and  decreases  slowly  with  distance  from  the  trunk  and 
is  shown  to  be  approximately  IdB  at  30  meters  from  the  trunk.  It  is  also  shown  that 
the  method  of  calculating  the  scattered  field  near  a  tree  trunk  at  normal  incidence 
can  be  used  for  oblique  incidence  of  finite  length  cylinders.  However,  its  applicability 
fails  for  large  angles  of  incidence  and/or  for  large  distances  from  the  cylinder  when 
a  direct  ray  between  the  sotirce  and  field  point  dominates  the  near-zone  scattered 
field. 

The  backscattered  field  of  a  tree  trunk  is  examined  next.  First  the  effect  of  the 
trunk-ground  interactions  (TGI)  is  examined.  The  modified  Geometrical  Theory 
of  Diffraction  is  used  to  find  the  scattered  field  for  this  geometry.  It  is  found  that 
the  scattered  field  for  horizontal  and  vertical  polarizations  is  approximatdy  7dB 
different  in  the  dihedral  region  between  the  trunk  and  the  ground.  This  is  explained 
by  the  decrease  in  the  reflection  coefficient  for  both  the  trunk  and  ground  dielectric 
surfaces  near  the  Brewster  angle  for  vertical  polarization.  The  Brewster  angle  for  the 
horizontal  polarization  isn’t  present  for  this  geometry  since  no  magnetic  materials 
are  used  in  the  modeling. 

To  determine  the  effect  of  the  scattering  from  branches,  the  dihedral  between 
the  brandi  and  trunk  is  examined.  The  purpose  is  to  find  a  levd  for  the  scattering 
from  a  branch  and,  therefore,  to  find  the  importance  of  branch  scattering  compared 
to  the  trunk-ground  interactions.  Therefore,  since  the  exact  level  of  the  scattering 
from  a  branch  isn’t  important  but  the  maximum  level  is,  a  worst  case  analysis  is 
examined.  This  allows  the  branch  geometry  to  be  greatly  simplified.  It  is  found  that 
the  scattered  field  in  the  dihedral  repon  of  the  branch  scattering  is  approximately 
20dB  lower  than  the  TGI  scattering  in  the  same  region.  The  maximum  level  of 
branch  scattering  occurs  for  normal  incidence  on  the  side  of  the  branch  and  is  less 
than  12dBsm.  But,  since  the  inddent  fidd  will  rarely  be  normally  inddent  on  the 
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branch,  the  branch  scattering  should  be  less.  The  branch  scattering  also  exhibits 
the  same  difference  in  the  RCS  for  the  dihedral  region  between  the  horizontal  and 
vertical  p<dazization8  as  the  TGI  geometry.  The  cause  for  the  difference  is  also  due 
to  the  change  in  the  reflection  coefficient  near  the  Brewster  angle  as  is  the  cause 
for  the  difference  in  the  TGI  geometry.  Since  the  scattering  from  a  single  branch  is 
much  less  than  the  scattering  from  the  TGI  interactions,  for  several  branches,  the 
level  of  branch  scattering  will  increase  but  the  TGI  will  remain  larger  for  a  single 
tree. 

In  condunon,  for  the  low  frequencies  of  interest,  simplificatioi  s  can  be  made  in 
modeling  the  scattering  from  the  trunks  and  branches  of  trees.  Also,  the  higher 
order  terms  between  branches  and  the  branch  and  the  ground  will  be  much  smaller 
than  the  main  branch  scattering  for  the  geometries  of  interest  and,  therefore,  are 
ignored.  The  leaves,  smaller  branches,  and  higher  order  terms  were  ignored  in  the 
analysis,  but  would  need  to  be  included  to  completely  model  the  scattering  from  a 
tree. 
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Appendix  A 

Reflection  and  Transmission  of 
plane  waves  in  Planar,  Stratifled 
(multilayered)  Medium 


DetailB  of  reflection  and  transmission  of  the  plane  wave  in  multilayered  medium  are 
presented  in  [10].  Here  a  summary  of  [10]  is  presented.  Suppose  a  plane  wave  in 
flee  space  has  oblique  inddence  on  a  plane  multilayer  consisting  of  N  homogeneous 
isotropic  slabs  as  shown  in  Figure  53.  Let  dn,  Mn,  and  Cn  represent  the  thickness, 
permeability,  and  permittivity  of  slab  n.  The  slabs  are  considered  to  have  infinite 
width  and  height  and  parallel  surfaces,  with  unbounded  flee  space  on  both  sides  of 
the  multilayer.  The  inddent  plane  wave  impinging  on  the  left-  hand  surface  of  the 
multilayer  is  pven,  in  the  T£  case  (i.e.,  perpendicular  polarisation)  by: 

(35) 


where  d  is  the  angle  of  inddence,  ho  =  2ir/A,  and  A  is  the  flee  space  wavdength. 
The  reflected  wave  is  given  by: 


(36) 


where  R  is  the  reflection  coeffident  of  the  multilayer.  The  transmitted  plane  wave 
on  the  right-hand  ride  of  the  multilayer  is  represented  by: 


(37) 
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whete  T  is  the  traasinission  coeffident  of  the  multil&yer.  The  field  in  each  layer  can 
be  regarded  as  an  infinite  series  of  the  plane  waves  bonndng  back  and  forth,  but  it 
is  more  convenient  (and  equally  valid)  to  consider  it  to  be  the  sum  of  only  two  plane 
waves,  one  travelling  outward  and  one  reflected.  In  layer  n,  for  example,  the  field  is 
represented  by: 

En  =  (A«e^-*  +  B„e-’-')  (38) 

Similaxly,  in  layer  n+1  the  electric  field  intensity  is  given  by: 

(39) 

The  boundary  between  layers  n  and  n  + 1  is  located  at: 


Xn  =  di  +  cfe  +  d3  +  ***  +  dn-  (40) 

By  enfordng  the  boundary  conditions  on  Eg  and  JTj,  at  r  =  Xn,  it  can  show  that: 

(41) 


■^+1 

Pn 

Qn' 

An' 

Bn+j 

Rn 

» 

. 

Bn 

where: 


2  \  PnTn-fl  / 
2  \  Pn7n+l/ 


(42) 

(43) 

m 


and: 


5„  =  i  fl  +  I  (45) 

2  \  fin'in+l/ 

The  propagation  constant  jn  for  layer  n  will  be  complex  if  the  medium  is  dissi¬ 
pative.  Both  the  real  and  the  imagery  parts  of  7n  will  be  positive.  If  layer  n  is 

a  losdess  medium,  7n  will  be  purely  imaynary.  The  wave  equation  is  employed  to 

obtain: 

7i3  =  iv/w’Pnei.  (46) 
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The  reflection  and  transmission  coeffident  of  the  multilayer  can  be  calculated  in 
a  systematic  manner  by  setting: 


i4o  =  1 


(«) 


and: 

flb  =  0  (48) 

and  then  by  using  the  recursion  equations  Equation  (41)  to  calculate  Ai,  Bu 
Bit  •••  >nd  Bn+i  in  that  order. 

From  Eqs.  (35)  through  (38), 


Eo  =  Aat+i, 

(49) 

_  Bjv+i 

(SO) 

Aat+i* 

(51) 

In  the  TE  case  the  constants  An  and  Bn  represent  the  electric  Add  intensities  of  the 
outgoing  and  reflected  wave  in  each  layer.  In  the  TM  case  (paralld  polarisation) 
the  solution  proceeds  in  the  same  manner.  The  equations  pven  above  apply  in  both 
cases  but  the  A„  and  Bn  represent  the  magnetic  Add  intensities  in  the  TM  case  and 
/tn+i  sad  Mn  must  be  replaced  with  Cn+i  and  Cn  in  Eqs.  (42)  through  (45). 

If  a  perfectly  conducting  sheet  is  placed  on  the  right-hand  surface  of  the  multi¬ 
layer  (i.e.,  on  the  z-y  plane),  the  solution  is  again  pven  by  the  equations  above  with 
the  exception  that  the  transmisnon  coeffident  T  is  not  cdculated  in  this  case,  and 
Eqs.  (47)  and  (48)  are  replaced  with: 


Ai  =  1  and  J9i  =  — 1  in  the  TE  case  (52) 

and: 

Ai  1  and  =  1  in  the  TM  case.  (53) 


Equatimu  (52)  and  (53)  axe  obtained  by  forcing  the  tangentid  electric  Add  intennty 
to  vanish  at  the  perfectly  conducting  plane. 


In  the  above  equations,  the  reflection  coefficient  R  is  defined  as  the  ratio  of  the 
reflected  wave  amplitude  to  the  incident  wave  amplitude  at  the  coordinate  ori^; 
that  is, 


n  g  (0,0.0) 
£•(0,0,0)’ 


for  the  TE  case. 


(54) 


and: 


flJ(0A0). 

“■'£•(0,0,0)’ 


for  the  TM  case. 


(55) 
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Appendix  B 

Scattering  of  Plane  Waves  by 
Stratified  (multilayered) 

Cylinders 

Details  of  the  scattering  of  phae  waves  by  moltilayeted  dielectric  cylinders  of  infinite 
length  with  nonnal  incidence  ate  presented  in  [10].  Here,  a  summary  of  [10]  is 
presented.  Consider  a  plane  wave  to  have  nonnal  incident  on  a  multilayered  dielectric 
qdiader  consisting  of  M  lossless  homogeneous  layers  as  riiown  in  Figure  54.  The 
qrfinder  axis  is  tahen  to  be  the  a  axis  in  the  rectangular  coordinate  system,  and  the 
z  axis  is  the  axis  of  propagation  of  the  incident  plane  wave.  The  permeability  and 
permittivity  of  layer  m  ace  denoted  by  pm  and  Cm,  and  the  phase  cmistant  is  pven 
by: 

km  =  V^/|^mtm  (56) 

In  the  TM  case,  the  dectric  fidd  has  onfy  a  s  component  ghren  in  layer  m  by: 

=  +  (57) 

nmO 

where  (p,^,s)  are  the  cylindrical  ooor^nates  and  Jn  and  rqnesent  the  Bessd 
and  Neumann  functions.  The  time  dependence,  e^,  is  implied.  In  the  TE  case  the 
sdution  is  obtained  from  the  equations  pven  here  by  interchanging  p  and  e  and  E 
and  Sf  where  S  represents  the  magnetic  fidd  intenrity. 
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Figure  54:  A  Mvhilajpered  Didectric  Cylinder  of  Infinite  Length. 
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TAe  coeffidentt  Amn  and  Bmn  must  be  determined  by  applying  the  boundary 
conditioins  of  E,  and  at  each  interface.  From.  Maxwell  equations,  the  ^  component 
of  the  magnetic  field  intensity  in  layer  m  is  pven  by: 

(58) 

naO 

where  J'  and  JV'  denote  the  derivatives  of  the  Bessel  and  Neumann  functions  with 
respect  to  the  arguments. 

Let  Rm  denote  the  outer  radius  of  layer  m.  Rrom  the  boundary  conditions  at 
the  inter&ce  between  the  layers  m  and  m+1,  it  is  found  that  the  wave  expansion 
ooeffidents  for  the  two  refpons  have  the  following  linear  rdations: 

“  ^mnAmn  4" 


where: 


^  (61) 
Kim  ~  lhnKn+lJn{,KiiRm)JniKm+lEm)]  (62) 


^mn  —  ■^Jjj^[/*«+l^'lVj^(^'Rm)«fi»(hm+l^m)  Mfi»^+l^»(^'®m)*fB(^+l^^)l 


If  the  coefficients  in  the  first  layer  (Ain  and  Bm)  were  known,  the  coeffidents  in 
tlm  remaining  layers  could  be  calculated  by  rang  £quadons(59)  and  (60),  tecurnvdy. 
To  permit  a  procedure  of  this  type,  define  a  set  of  normaHaed  ooeffidents  and 


the  proportionality  constant  as  follows: 


rr'  ' 


Amn  —  ‘^n^fitn 


(65) 


and: 


Bmn  —  ■K^B-Brnn* 


(66) 


M^th  no  loss  of  generality,  Let: 

Ain  =  1 


(67) 


U  the  center  layer  is  a  dielectric  medium,  the  field  mnst  be  finite  at  the  origin  and: 


=  0. 


(68) 


The  normalised  coeffidents  also  obqr  the  same  recursion  formulas  and  we  can  now 
calculate  and  m  that  order. 

The  field  in  the  exterior  free-space  repon  is  pven  l^: 


=  f:i(-jr«./.(*o/>)  +  (69) 

MaO 


where  H^\hop)  represents  the  Hankel  function  and  ko  is  the  phase  constant  of  free 
space.  The  first  series  in  £qnation(69)  rq>resents  the  inddent  plane-wave  fidd,  and 
the  second  series  is  the  scattered  field  which  contains  outward-travding  waves  only. 
The  function  Cn  is  unity  if  n=0,  and  en=2  if  n  is  greater  than  sero. 

The  scattering  coeffidents  for  the  exterior  re|pon  are  pven  by: 


—  3Bt4^\^  — 


^S/+l.n  jA,n+l,n 


(70) 


This  otnnpletes  the  solution.  Equation(57)  can  be  used  to  cdculate  the  fidd  at  any 
point  in  the  dielectric  ^rlinder,  and  £quations(69)  and  (70)  are  used  to  cdculate  the 
field  at  ai^  point  outside  the  cylinder. 

In  the  ease  of  a  perfect^  ocmducting  cylinder  with  one  or  more  dielectric  layers, 
let  a  be  the  radius  of  the  conducting  cylinder  and  Ri  the  outer  radius  of  the  first 
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dielectric  layer.  The  above  equations  again  give  the  solution  if  £quation(66)  is 
replaced  with: 


for  the  TM  case  and: 


B[n  = 


B[n  = 


Jnjkia) 

>«(ifeia) 


(71) 


(72) 


Kikia) 

for  the  TE  case.  Equations(71)  and  (72)  are  obtained  by  setting  £,=0  or  at 
the  conducting  surfisce. 

When  a  plane  wave  is  incident  on  a  scattered  cylindrical  structure  of  infinite 
length,  the  distant  scattering  pattern  is  described  by  the  echo  width  which  is  defined 


w 


(73) 


tr2D  —  2tp 

The  incident  electric  fidd  E*  is  taken  to  have  unit  magnitude.  The  scattered 
dectric  fidd  for  large  distances  from  the  cylinder, 


where: 


(74) 

-B'm+i^ 

(73) 

Substituting  Equations(74)  and  (75)  into  Equation(73)  the  Instatic  echo  width 
of  the  multilayered  dielectric  cylinder  is  (pven  by: 

|2 


2Ai~ 
ff2D  =  - 

W 


^enI>nCOBn^j 
fisO 


(76) 


Appendix  C 

Scattering  &om  Multijoined 
Cylinders 


This  appendix  snmmaxises  the  solution  for  the  mnltyoined  cylinder  geometry  solved 
by  N.  Akhter  [25].  The  geometry  for  the  mnltijmned  cylinders  is  shown  in  Figure  55. 
The  Uniform  Theory  of  Difiaction  is  used  to  calculate  the  £sr*aone  scattered  field 
from  this  geometry.  The  solution  in  based  on  the  modification  of  a  2*D  dihedral 
solution  using  UTD  by  N.  Alchter  [25].  The  solution  is  modified  for  dielectric  material 
using  the  dielectric  wedge  diffraction  coefficient 
in  Equation  (6). 

The  interactions  needed  to  calculate  the  backscattered  field  in  the  UTD  solution 
are  shown  in  Figure  56.  The  diffracted  field  for  each  interaction  term  is  calculated 
similar  to  the  interactions  for  the  trunk-ground  interaction  geometry  in  Chapter  6. 
For  the  angle  difficaction  D  terms,  the  calculation  is  the  same  as  in  Chapter  6  except 
now  there  axe  2  cyfinders.  The  spread  factors  for  the  D  terms  shown  in  Equation  (2) 
axe  calculated  nmilady  to  the  D  terms  in  Chapter  6  using  Equation  (5). 

The  diffracted  fidd  for  the  RD  and  DR  terms  b  calculated  iriTniln.T  to  the  RD  and 
DR  terms  in  Chapter  6.  The  difference  b  that  the  multyoined  cylinder  geometry 
consists  of  2  cylinders  versus  1  cylinder  for  the  trunk-ground  interaction  geometry. 
Thb  diffexenoe  allows  rays  to  reflect  and  diffract  from  both  cylinders.  Therefore, 
terms  induding  reflections  from  cylinder  ^2  and  diffractions  from  cylinder  need 
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to  be  calculated  in  addition  to  the  tenns  including  reflections  from  cylinder  #1  and 
diftactions  flrom  cylinder  #2.  The  difl^action  coefficients  for  the  RD  and  DR  fields 
are  similar  to  the  diifraction  coefficients  in  Chapter  6.  However,  since  both  surfaces 
are  now  curved,  the  spreading  factor  will  be  different  than  the  spreading  &ctor  used 


in  Chapter  6. 

The  spread  factors  for  the  RD  and  DR  terms  shown  in  Equation  (2)  are  calculated 
using  Equation  (5).  For  the  RD  and  DR  terms,  the  spread  factors  will  be  the  same 
since  the  rays  travel  along  the  same  ray  path.  The  spread  &ctor  for  the  ray  diffracted 
from  cylinder  #1  and  reflected  from  the  nde  of  cylinder  #2,  the  spreading  factor  is: 


where: 


o  +  (a  +  p")2sin(a  —  fi) 


A  ^  » 

^  sin  d  +  8in(2a  —  B)' 


In  the  above  equations,  a  is  the  (^linder  radius  and  a  is  the  dihedral  an^e  between 
the  cylinders.  The  spread  factor  for  the  ray  dii&acted  from  cylinder  #2  and  reflected 
from  the  nde  of  cylinder  #1,  the  spreadii^  factor  is: 


a  +  (a  +  p“)2rin^ 


sin(a  -f-  ^)  +  sin(a  —  B) 

The  final  interaction  of  interest  is  the  Difl^acted-Reflected>Diffracted  (DRD)  term 
which  is  calculated  similar  to  the  previous  interaction  terms.  The  spread  factor  for 
the  ray  diffracted  from  cylinder  #2,  reflected  from  blinder  #1,  and  then  diffracted 
from  cylinder  #2  again  is  given  by: 


JhiiiDa  _  [  _ Pifi  d 

l»W  +  ')  W  +  •Xpi  + 


(81) 
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